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ABSTRACT 

The problem of chemo-photometric evolution of late-type galaxies is dealt with relying on 
prime physical arguments of energetic self-consistency between chemical enhancement of 
galaxy mass, through nuclear processing inside stars, and luminosity evolution of the system. 
Our analysis makes use of the Buzzoni template galaxy models along the Hubble morpholog- 
ical sequence. The contribution of Type 11 and la SNe is also accounted for in our scenario. 
Chemical enhancement is assessed in terms of the so-called "yield metallicity" (Z), that is 
the metal abundance of processed mass inside stars, as constrained by the galaxy photometric 
history. For a Salpeter IMF, Z cx <°'^'^ being nearly insensitive to the galaxy star formation 
history. The ISM metallicity can be set in terms of Z, and just modulated by the gas fraction 
and the net fraction of returned stellar mass (/). For the latter, a safe upper limit can be placed, 
such as / < 0.3 at any age. 

The comparison with the observed age-metallicity relation allows us to to set a firm upper 
limit to the Galaxy birthrate, such as 6 < 0.5, and to the chemical enrichment ratio AY/ AZ < 
5. About four out of five stars in the solar vicinity are found to match the expected Z figure 
within a factor of two, a feature that leads us to conclude that star formation in the Galaxy must 
have proceeded, all the time, in a highly contaminated environment where returned stellar 
mass is in fact the prevailing component to gas density. 

The possible implication of the Milky Way scenario for the more general picture of late- 
type galaxy evolution is dicussed moving from three relevant relationships, as suggested by 
the observations. Namely, i) the down-sizing mechanism appears to govern star formation in 
the local Universe; ii) the "delayed" star formation among low-mass galaxies, as implied by 
the inverse 6-Afgai dependence, naturally leads to a more copious gas fraction when moving 
from giant to dwarf galaxies; Hi) although lower-mass galaxies tend more likely to take the 
look of later-type spirals, it is mass, not morphology, that drives galaxy chemical properties. 
Facing the relatively flat trend of Z vs. galaxy type, the increasingly poorer gas metallicity, as 
traced by the [O/H] abundance of Hll regions along the Sa Im Hubble sequence, seems to 
be mainly the result of the softening process, that dilute enriched stellar mass within a larger 
fraction of residual gas. 

The problem of the residual lifetime for spiral galaxies as active star-forming systems 
has been investigated. If returned mass is left as the main (or unique) gas supplier to the 
ISM, as implied by the Roberts timescale, then star formation might continue only at a max- 
imum birthrate bmax ^ //(I ~ /) ^ 0.45, for a Salpeter IMF. As a result, only massive 
(A^/gai ^ 10^^ -^-^o) Sa/Sb spirals may have some chance to survive ^ 30% or more beyond 
a Hubble time. Things may be worse, on the contrary, for dwarf systems, that seem currently 
on the verge of ceasing their star formation activity unless to drastically reduce their apparent 
birthrate below the 5,nax threshold. 

Key words: Galaxy: disc - Galaxy: evolution - galaxies: abundances - galaxies: ISM - 
galaxies: spiral 



1 INTRODUCTION 

A fair estimate of "galaxy metallicity" has been widely recognized 
as a basic issue when assessing the distinctive evolutionary prop- 



erties of late-type galaxies (LTGs), including all disk-dominated 
systems along the Hubble morphological sequence. However, since 
both spiral and irregular galaxies consist of a composite stellar pop- 
ulation, the question cannot easily be settled, as it depends on the 
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specific criterion we want to adopt to identify a value of Z (or 
equivalently of [Fe/H]) ideally representative of the stellar system 
as a whole. 

For example, the youngest stars closely track the chemical 
composition of gas clouds within a galaxy, and we could therefore 
refer to their value of Z as the galaxy actual metallicity [Zo)- On 
the other hand, this value may not account for the composition of 
earlier generations of stars, that however still provide a sizable con- 
tribution to galaxy total luminosity. In this regard, one would better 
like to consider an effective metallicity estimate (Zeff ), by weight- 
ing the composing stellar populations with their actual photometric 
contribution or with the number of star members in each group. 

T hese arguments y et emerge "in nuce " from the semi nal anal- 
ysis of lSchmidtl (l963), but onlv lArimoto & Yoshiil ( Il987i) exten- 
sively dealt with the problem of a fair estimate of the effective 
metallicity in their LTG models using bolometric luminosity as 
a weighting factor to account for the relative contribution of the 
composing stellar populations. With this criterion, by considering 
a collection of theoretical si mple stel l ar pop ulations (SSPs) with 
a spread in the values of Z, iBuzzonil ( 1 19951) showed that the in- 
tegrated spectral energy distribution (SED) for the aggregate still 
resembled the spectrum of jus t one SSP with log Z^a ~ (log Z). 
From her side, iGreggiol (Il997h further refined these arguments, at 
least in the early-type galaxy framework, emphasizing the differ- 
ence between average methods such as to lead to [< Fe/H >] vs. 
< \Fe/H] >. A direct estim ate of both Zo and ZcS is provided by 
lArimoto & Jablonkal ( Il99ll, hereafter AJ91) in their LTG models. 
As expected, the effective metallicity is found to be systematically 
lower by a 20-30% with respect to the actual one. 

Apart from any operational definition, however, one has to 
remark that the classical approach to galaxy chemical evolution, 
usually pursued in the literature, mainly proceeds from a dynami- 
cal point of view, where chemical composition is tightly related to 
the kinematical properties of the different stellar populations shar- 
ing the gala xy morphological design. Sandage & Fonts ( 1987) and 
lKoeDpen & Arimotd (l990|) are striking examples, in this regard, 
of the empirical and theoret ical focus to t he problem, respectively. 
After all, just on these basis lBaad3 ( 1 19441) led to his first definition 
of Pop I and II stars. 

However, an even more direct physical argument ties metal 
abundance to galaxy luminosity evolution. For any gram of met- 
als eventually produced in a galaxy, in fact, one gram of Hydro- 
gen has been burnt at some time releasing a fixed amount of en- 
ergy. As a consequence,a straightforward relationship must be in 
place between the actual amount of enriched mass in a galaxy 
(Myz), in the form of metals plus Helium in excess to the primor- 
dial abundance, and the past luminosity evolution. The importance 
of this constraint, often overlooked even in the recent literature, 
has been first raised by WirthI ( Il98ll) in the contex t of ea r ly-typ e 
galaxy evolution, and more explicitely assessed bv lPagell ( Il99% . 
This argument also found a more definite theoretical settlement 
through the so-called "Fuel consumption theorem", as proposed by 
[Renzini & Buzzoni (1983, 1986). 

Given the range of interesting implications for the study of 
galaxy evolution, we want here to further investigate the problem of 
the energetic constraint to chemo-photometric evolution of galax- 
ies, with special attention to the more composite scenario of late- 
type systems . Our ana l ysis \yill rely on the template galaxy models 
developed in lBuzzonil jlOOSl 12003) . To better catch the real inten- 
tion of this paper, we would like to warn the reader that our dis- 
cussion will proceed as far as possible on analytical bases. This 
will forcedly require some educated simplifications to the problem 



(after careful physical assessing however of the neglected compli- 
cations), with the aim of more neatly let to emerge the leading 
relationships that govern galaxy evolution, and in particular the 
way these systems (and the Universe as a whole) constrained their 
chemical composition at the different cosmic epochs. 



2 THEORETICAL FUNDAMENTALS 

In order to better single out the relevant physical boundaries in- 
volved in our problem, we need to consider basically the two main 
photometric contributions, that constrain overall galaxy luminosity 
evolution. From one hand, in fact, we can recognize the "smooth" 
photometric contribution by stars along their quiescent (namely 
"hydrostatic") evolutionary stages, according to a given IMF. On 
the other hand, one has also to consider the Supernova bursting 
events, that interleave galaxy evolution with short but extremely 
relevant episodes of mass processing. 

Following Clayton (1983), H-buming cycles release 26.731 
MeV per *He particle, while a supplementary amount of 7.274 
MeV per ^^C nucleus is produced by the 3q He-burning cycle. 
Overall, the energy budget of the nuclear processes that lead to Car- 
bon synthesis in stars reaches a total of 

W^, = 6.98 10^* erg g"\ (1) 

This specific energy output accounts for the "quiescent" luminosity 
emission of a galaxy provided by the non-explosive evolutionary 
stages of stars along the entire IMF mass range. 

A supplementary important source of energy and ISM metal 
enrichment should however comprise also the SN explosive events. 
A quite different elemental pattern characterizes the ISM chemi- 
cal enrichment in case of SNe la, meant to be the main synthe sizers 
of the Fe-Ni elemental group dChevaliedl 19761; iNomotolll 980 ). and 
SNell, which stem on the contrary from the core collapse of high- 
mas s stars and mainly supply Oxygen and other heavier a elements 
(e.g. lMatteucci & Greggi(jl98^jTrimblell99"ll) . Nevertheless, both 
these events may actually be treated coherently in terms of specific 
energy release as they take in charge the post-Carbon nucle ar burn- 
ing up to Iron synthesis (and beyond). Again, following iClavtonI 
( Il983[) . whatever the chemical path of nuclear reactions, the atomic 
binding energy difference implies to release 60.605 MeV per ^^Fe 
nucleus or, equivalently, 

W^SN = 1.03 10'%rgg-\ (2) 

Compared to "quiescent" nuclear reaction chains, for fixed 
amount of released energy, explosive nucleosynthesis is therefore 
roughly a factor of seven more efficient in the whole Z enrichment 
of the Universe, selec tively involving, however, only the heaviest 
metal supply (see, e.g. iMatteuccill 2004 for a brief but timely focus 
on this subject). 

2.1 Processed mass by "quiescent" luminosity evolution 

Let us explore first the general case of a system consisting of a com- 
posite stellar population (CSP), as a result of a roughly continual 
feeding of fresh stellar generations along the whole photometric 
history. The nature itself of the problem suggests to consider the 
total bolometric luminosity of the CSP as a convolution integral of 
several simple stellar populations (SSPs) continually distributed in 
time according to a given star formation rate (SFR), namely 

£csp = /^ssp(r)® SFR dr. (3) 
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In iBuzzonil j2005l) we have shown that SSP bolometric lu- 
minosity (^ssp) smoothly evolves with time according to a sim- 
ple power law such as Issvit) oc t~°'\ for a Salpeter IMF, the 
power index a = 0.77, is virtually independent from the SSP 
chemical composition^ Our template galaxy models considered a 
disk SFR wi th a characteristic birt hrate, b = SFRo/ (SFR) (see, 
for instance Miller & Scald 1979|), depending on the morpholog- 
ical type l lKennicutt. Tamblvn, & CongdoiJI 19941 ). The distinctive 
gas — > star conversion efficiency implies for the star formartion 
rate a power-law time dependen ce such as SFR(t) « t~'^, so that 
6 = (1 - /3) (see lBuzzonilbOOl for a full discussion). 

Within this framework, the total luminosity of a CSP at time t 
becomes therefore 



£csp(t) = C 



- it- 



dr, 



(4) 



with a normalization constant C depending on the total mass of 
the aggregate. Providing that both a < 1 and /3 < 1, the integral 
equation has the Euler Beta function as a straightforward analytical 
solution. More explicitely, it can also be written as 

r(i-Q)r(i-/3) 



£csp(t) 



r(2 - Q - p) 

being F the Gamma function and 7 - 
ence epoch to, one can simply write 



Ct'-'', (5) 
Q + /3. Compared to a refer- 



1-7 



(6) 



If we integrate the CSP emitting power jCcsp over the entire 
galaxy life, until time T, this eventually gets the global energy £ 
released by the nuclear processes to sustain galaxy luminosity until 
that age. 



SiT)^ f j!:csp{t)dt = £o 

Jo V'-o 



(7) 



where, £0 = [Ccsp{to) x to]/ (2 - 7). 

Clearly, at any time the latter quantity must balance the 
amount of processed mass Myz. With respect to the total CSP stel- 
lar mass (that is the total mass ever converted in stars, M* (T) — 
SFR(t) dt), we have that 



M* 



Wr^£{T) 
Jg SFR(f) dt 



£0 



W,Mo 



(8) 



where W, is the energy-mass conversion factor set by the nu- 
clear burning, according to our previous argumentsjf] Again, we 
set Mo = M*{to) = Jg*° SFR(r) dr for notation convenience. 
Interestingly enough, eq. ((8) shows that the way metals enhance 
with time in a CSP does not depend on the SFR details, being 
fully modulated on the contrary by the IMF slope, which even- 
tually constrains the CSP luminosity evolution through the power 
index a. Furthermore, in the same equation, note that a parameter 
J- = {£o/Mo) /W, can also be regarded as a "burning efficiency 
factor", which is a measure of how deeply processed mass has been 
exploited to produce luminosity. 



Our SSP evolution assumes a ISalpeteJ fT95^ IMF such as M{M) oc 
M"", with the power index s = 2.35 along the accounted stellar mass 



range from 0. 1 to 120 . See lSuzzonil j2005il2002l) for further details. 
^ Note that W, directly relates to the total efficiency e of nuclear reac- 
tions. According to the Einstein equation, i.e. E = e AMc^, we have that 
Wt = £/ A.M = ec^. This eventually leads to a mass — ^ energy conver- 
sion efficiency of e = 0.0078. 



Table 1. Reference values at t = 15 Gyr for a lO^lM, 
with changing star-foiTnation properties'"' 
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0.110 
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0.108 
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1.02 


0.98 


7.33 


0.105 


-0.2 


1.2 


0.89 


1.12 


7.21 


0.103 


-0.4 


14 


0.80 


1.25 


7.06 


0.101 


-0.6 


1.6 


0.72 


1.38 


6.94 


0.099 


-0.8 


1.8 


0.66 


1.51 


6.84 


0.098 


-1.0 


2.0 


0.61 


1.64 


611 


0.097 


Assuming SFR 


oc t ^ , and a birthrate 6—1 


-f) 





For a 1 Mm CSP 



Burning efficiency factor T = (So/Mo) K ^ 

2.1.1 Trading mass vs. luminosity: the galaxy M/L ratio 

Reference stellar mass (Mo) and the corresponding CSP luminosity 
(>Ccsp (to)) are, of course, strictly related depending on the SFR de- 
tails. Relying on the Buzzoni (2005, 2002) theoretical framework, 
we computed in Table [T] the expected value of the M/L ratio at 
15 Gyr for different CSPs (col. 3 in the table) along a range of 
values of the SFR power index, /3 (col. 1) and stellar birthrate b 
(col. 2). These values are representative of the full range of scenar- 
ios that characterize disk evolution for the different LTG models. 
From the values of M/L one can easily derive the expected CSP 
bolometric luminosity, £csp(io), as well as the value of £o at the 
reference age (col. 4 and 5, respectively, of Table [T). Finally, the 
burning efficiency factor, F is also computed in column 6. 

As both the total stellar mass and luminosity respond in the 
same way to the SFR, one has that 



M 



M_ 



(9) 



that is time evolution of galaxy M/L ratio is a nearly universal law, 
that only depends on the observing wavelength and, for bolometric 
and a Salpeter IMF, scales aproximately as M/L oc t'^''^. 

The SFR enters, on the contrary, by setting the absolute value 
of M/L, at a given age; for example, when compared to a simple 
star burst, a smooth star formation will result, on average, in a larger 
number of young bright stars of high mass at a given age. This 
means, in general, that M/L decreases with decreasing /3. In this 
regard, a convenient fit to the data of Table [T|provides: 



(M/L) I 



-f 0.19, 



(10) 



(0.23 + b) 

with an rms of ±0.04 on the predicted ratio. 
2.2 Supernova processing 

The "quiescent" evolution of galaxy luminosity, as described by 
previous equations, does not take into account the extra-energy re- 
lease (and the corresponding extra-metal production) by the SN 
events. As well known, supernovae display a composite genesis. 
From one hand, in fact. Type II SNe basically deal with the C/0 
core collapse of individual stars with conveniently high mass such 
as to exceed the Chandrasekhar limit near the end of their photo- 
metric evolution. On the other hand, a much more composite family 
of objects (shortly grouped into Type la SNe) has to be reconducted 
to binary-star evolution. In this case, the key mechanism deals with 
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the mass transfer between the system members, where a primary 
(originally more massive) star, in its way to the final white-dwarf 
(WD) death, acts as a fresh-matter supplier to boost evolution of 
the secondary member. The latter star will therefore be prompted, 
at some stage, to feed back the WD with fresh fuel until triggering 
(under appropriate physical conditions) the fatal detonation on the 
WD surface. 

As well known, the complexity of this scenario, characterized 
by a number of unknown or poorly constrained physical parame- 
ters (i.e. orbit parameters, total mass of the system, and mass ratio 
of the members), is the main source of uncertainty even in more ex- 
plicit and sophi sticated models of galaxy chemical evolution (see, 
in particular, Matte ucci & Greggioll 19861 for a brilliant assessment 
of the problem). All the more reason, let us try here a simplified ap- 
proach, just to single out the few leading quantities that constrain 
our relevant physical output. 



2.2.1 TypellSNe 

Considering first the core-collapse events, origi nating by the explo- 
sion of stars with A/up > 8 ± 1 M© (see, e.g. ISmartt et J]|2009l) 
one has that the expected total number of events per unit SSP total 
mass can be written as 
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Figure 1. The expected fraction of enriched stellai' mass provided by the 
SN II events according to three different IMF power-law indices (Salpeter 
case is for s = 2.35), and for a different value of the SN triggering mass, 
A/up. The consensus case is mai'ked by the big dot on the Salpeter curve. 
Vertical en'or bais derive from eq. (1141 . 



As a consequence, the appearence of SNiI on the scene of galaxy 
evolution just enters with a straight offset to metallicity|f| 



nsNii(SSP) = 



2 - s 



[M^-]jJu°p 



1120 

lo.i 



(11) 



where we assume a power-law IMF according to footnote [T] re- 
marks. For a Salpeter case (i.e. s = 2.35), and in force of the fact 
that galaxy SFR does not change much on a time-scale comparable 
with the lifetime of st ars with M ^ A/up (namely a few 10*^ yr at 
most, see e.g. eq. 3 in lBuzzoriill2002h . we have 



nsNii(t) = 0.126 A/u"p''^^ SFR(t) 



(12) 



For A/up = 8A/q, this leads to a rate of 7.6 10"^ SNiI events 
per year for a SFR — 1 Mq yr~^. This figure perfectly matches 
the empirical rate of (7. 5 ± 2. 5) 10~^ independently obtained by 
IScannapieco & BildstenI ( |2005|) on the basis o f the observed SN 
rate density within z ^ 1 bvl Dahlen et aL U2004t) compar ed against 
the cosmic SFR density as measured by Giavalisco et al. ||2004). 



If ES 



is the mean released energy associated to each SN 



event (in unit of 10'' erg), and recalling the SN specific energy 
release as previously discussed, then the fraction of enriched stellar 
mass (Myz) along time results 



M* 



10^ 



WsnMq 



: 0.061 A/u"p'-^^ E^l"', 



(13) 



being Wsn the energy-mass conversion factor set by the explo- 
sive nuclear burning as in eq. |2] and Mq the mass of the Sun. 
The output of eq. ( I13t is displayed in graphical form in Fig. [T] 
exploring also differen t IMF slopes. For A/up = 8 ± 1 Mo, and 



Ejl"' = 1.3 ± 0.3 (Thielemann, Nomoto. & Hashimotdll99d : 
lHashimotdri995l :l Kasen & Plewa.2007. ) one derives 



M- 



SNII 



M* 



= 0.0047^2?- 



(14) 



As for eq. ((D, note that both Myz" and M* depend on / SFR{t)dt, 
so that the Myz''/M* ratio does neither depend on time nor on the 
SFR details, being fully constrained by the IMF slope alone (and 
by the preferred value of A/up, as well), as displayed in Fig. [T] 



2.2.2 Type la SNe 

From the theoretical side, i Greggio & Renziiiil (Il983l) have first ar- 
gued that the SNia population originates from the mix of two stellar 
components. In fact, the SN detonation may set on both on He-rich 
WDs (mainly coming from the evolution of relatively low-mass 
stars, about 2-3 Mq), and on the surface of C/O-rich WDs (the final 
output of more massive stars, about 5-8 Mej^This difference, to- 
gether with the binary nature of the physical mechanism that drives 
the explosion, makes the total SNia rate to depend (at least) on two 
reference parameters. From one hand, simple evolutionary argu- 
ments on the stellar lifetime sets a minimum characteristic "delay 
time" (AsNi) for the SN burst, which may follow the star forma- 
tion ev ent (and the first SNiI burst) by a few 10^ yrs for CO-rich 
WDs ( Tomambe & Matteucci 1986) up to several 10** yrs for He- 
rich WDs (iGreggio & Renzinl.I983i) . On the other hand, the broad 
range of geometric and physical parameters of the binary system 
may even further displace in time the SN event, blurring the con- 
tribution of the two feeding channels into a composite delay-time 
distribution, namely 0(a) t hat we eventually have to deal with. 

From the observing side, iMannucci et al.l fcoOSh have firmly 
demonstrated that SNia rate among spiral galaxies tightly re- 
lates to the corresponding rate of core-collapsed SNe, the lat- 
ter meant to closely trace the galaxy SFR. This feature evidently 
sets im portant contrains to the Asni distribution, leading iGreggid 
l l2005h to conclude that, in a SSP "the distribution function of 



^ Given the special nature of the explosive SN nuclesynthesis, the reported 
j^SNll value should account for the emitted luminosit y, as obtained for 
instan ce from the time integration of the SN lightcurve jMilone & Milon j 
Il988l e.g.), plus the kinetic energy associated to the SN burst. The latter 
contribution is actually by fai' the prevailing one, and it exceeds by about 
two o rders of magnitude the emitted luminosity itself jWooslev & Jankal 
|2005|) . 

^ In any case, the upper mass of SNi progenitors must be edged by the 
on-set mass (Mup) of core-collapsed SNell. 
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the delay times decreases with time. This means that the major- 
ity of SNia precursors are relat ively short lived". On this line, 
iMaoz. Sharon.& Gal-Yaii] ( 1201 O h argued that over 50-85% of the 
whole SNia population in a SSP are expected to explode within the 
first Gyr from the starburst. On the other hand, the residual pres- 
ence of Type la SNe even among quiescent early-type galaxies, 
that is in lack of any evident star formation activity, calls for a few 
very delayed events and therefore for a skewed lo ng-term tail of 
the delay-time distribution ( Dallaporta 1973; Mannucci et al .120051 : 



Im Sd Sc Sb Sa 



IPanagia. Delia Valle. & Mannuccill2007l : iBrandt et aljboiol) . The- 
oretical fig ures for th e SSP case indicate that (jy^i^) t~^, with 
?7 ^ 1.2 jMaoz. Sha ron.& Gal- Yam 20ig), or ~ 2.3, as in the 



iGreggiolJioOSh or 
Similarly 



Matteucci & Recchil i20(m models □ 
to what already attempted 



. . by 

IScannapieco & BildstenI 1 2005h . we adopt for our discussion 
the iMannucci et alj ( |2005[) empirical parameterization for the 
observed SNia rate (see eq. 2 therein). One can therefore envisage 
a prevailing term directly r elated to galaxy SFR and an "ex tended" 
SN fraction (to recall the IScannapieco & BildstenI |2005| original 
notation), that collects the most delayed events and can therefore 
be meant to scale at any time with J SFR(t)dt, that is with the 
total stellar mass of a galaxy. 

For the first component, the direct SFR dependence is obtained 
through the core-collap sed SN rate as from eq. (1121 assuming, with 
iMannucci et al.l j2005l) , that 



0.008 



0.006 



|> 0.004 



0.002 
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15 Gyr 
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> 
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Figure 2. The chemical enrichment from SN contribution according to the 
different components. Note the prevailing role of core-collapsed objects (i.e. 
SN II) along the early stages of galaxy evolution (t < 50 Myr, see labels 
to the right axis). The total SN la component adds then a further 0.2 dex to 
galaxy metallicity, mostly within the first Gyr of life. 



nsNi 



0.35 nsNii- 

± 8 



By relying then on eq. ( |13t , this leads to 



M* 



0.35 



ipSNII 
^51 



Mi 



(15) 



(16) 



where the associated energy for a typical SNi event (again, 
in unit of lO^'^ erg) is of the same order than for SNell, 



name ly Eif 



1984 



1.5 ± 0.5 dNomoto. Thielemann. & Yokoil 



Khokhlov, Mueller, & Hoeflich| 

2002, 



19931: 



Reinecke. Hillebrandt. & Niemeyer 
2003h . 



iHillebrandt et all 



However, the l.h. term of eq. ( 116) has to be regarded as a some- 
what asymptotic value for the enriched mass fraction as it considers 
virtually no delay in the SNia appearence, being both enriched and 
total stellar mass perfectly tuned with SFR. If we set SNela to come 
on the galaxy stage just after the first SNiI metal "glitch", say for 
AsNi ^ 50 Myr or so, i.e. the reference lifetime for a 5 Mq star, 
then one has to scale previous equation by a factor (1 — r''"^), with 
T — t/AsNi. After this correction, we eventually obtain 



M- 



SNI 
YZ 



M* 



0.0019 



(17) 



where the error bars have been conservatively computed by loga- 
rithmic error propagation. 

To this contribution one has further to add the long-term "ex- 
tended" SN compone nt whose number of events, according to 
IMannucci et al.l j2005l) . amounts to 



^ Note that, in all cases, rj > f). This assures 0(a) f^de with time at a 
quicker rate compared to the reference SFR for our model galaxies. In force 
of this argument, we can therefore expect the bulk of SNia to safely trace 
the SFR all the time. 



nsNi = 4.7 10- 
± 9 



(18) 



This rate is nicely confirmed, within the errors, also by 
ISullivanetai] ( l2006l) . 

In case of a star-forming CSP with power-law SFR, the mass- 
specific event number that occurred up to age T can be written as 



M* 



cxt 

= "-SNI 



/I/oSFR(r)drdt 
XrSFR(r)dr 



A' ext. 
SNI ?^SNI 



rfl 



(19) 



2-/3 ■ 

The corresponding amount of processed mass directly derives 



M* 



iVs 



M* 



T-,SNI 10 



WMq 



By replacing the relevant quantities, we eventually obtain 



M5 



M* 



AsNI T [1 



-(1+6)1 



2.9101^ 



l + h 



(20) 



(21) 



providing to express the delay time Asni in years. 

In conclusion, the relative amount of processed mass related 
to the whole SN activity in a CSP of SFR power-law index /3 can 
be summarized as 



lYZ 



M* 
with 



0.0047 + 







0.0019 yl(T, Asm, &) 



A(r,AsNi,fo) = (l-r-'')+ ^'^^ 



for r < 1 



for T > 1, 



-(1+6)1 



(22) 



(23) 



5.5 10"' l + h 
In Fig.[2|we displayed the relevant output in graphical form, 
for different ages and values of the SFR power index (3. In the plot 
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we also marked the representative values of /3 for spiral g alaxies of 
differe nt Hubble morphological type, as derived from the iBuzzonil 
l l2005b templates. 



2.3 Metals versus Helium: the enrichment ratio 

With the help of data in Table [T] through eq. (|6) and ([8}, it is im- 
mediate to derive both disk luminosity (i.e. £csp) and the amount 
of chemically enriched mass (Myz) along galaxy life. 

One may even rely on the expected enrichment ratio 
R = l^Y/AZ to single out the fraction of Helium and 
heavy elements eventually produced. The R parameter can 
be tuned up either empirically, for example vi a the observa- 
tion of Galactic and extragala ctic Hll clo uds fcegueux e t al.l 
1979; Seii-ano & Peimbert 198j 
Peimbert. Peimbert. & Luridiana 



Paget etal.lll 992: Maciel 2001j; 



20021; ICasagrande et al.l 1200 



or theoretically, from stellar evo lution theory dMaeded 11993 ; 
iMarigo, Chiosi. & Kudritzkill2003h . Some unsettled discrepancies 
still exist, in this regard, with theory easing in general a more copius 
production of metals within the stellar population, thus implying a 
systematically lower value of R compared to the the observational 
evidence. As first pointed out by ,Mall ik & Mallik ( 1985), this dis- 
crepancy might be reconciled from the theoretical side by acting on 
the IMF upper edge, such as to decrease the expected contribution 
of the high-mass stellar component. AUover, an admitted range be- 
tween R = 3±1 is accepte d in the literature (e.g. Tosi 1996; M aciell 
l200ll ; lFukugita & Kawasaki 2006 ). and will be considered here. 

We could therefore introduce the concept of galaxy yield 
metallicity and Helium content defined, respectively, as 



Z = 



(A/* , 



M, 



(24) 



y = Yp + RZ, 



where y — Yp is the Helium abundance, in excess to the primordial 
content (Yp). 



3 RELEVANT MODEL OUTPUT 

The "yield" abundance of eq. ( 124b can be regarded as the actual 
mean composition of the bright stars in a CSP at a given epoch. Re- 
alistically, [Z, y) give a measure of the maximum amount of ISM 
contamination in the evidently extreme case of "diluting" the entire 
stellar mass into a vanishing residual fraction of fresh primordial 
gas. Accordingly, Z might also be regarded as the maximum value 
allowed to Zo, for newly bom stars. 

Clearly, the real process that leads to ISM chemical enrich- 
ment follows from a much more composite scenario, being the re- 
sult of a continual and subtly e ntangled interplay between gas and 
stars inside a galaxy ( lTinslevl[l98a) . Actually, this is the central 
subject of the many important theoretical contributions to the study 
of chemical evolution of galaxies that have been succeeding along 
the last few decades (see, e.g. Matteucci 2003, for an exaustive in- 
troduction to the subject). 

Nonetheless, even facing such a sophisticated reference 
framework, our simplified approach to chemical enrichment may 
still be worth of some attention, as definition of yield metallicity 
provides us with a straightforward and physically inherent tool to 
assess the maximum metal abundance that can be reached inside 
a galaxy at any time. With the study of [Z, y) we therefore place 
a firm constrain to galaxy chemical history in terms of the experi- 
enced luminosity evolution. A number of interesting considerations 



Table 2. Photometric properties, reference M/ L ratio and yield metallicity 
at t = 15 Gyr for galaxies of different morphological type'"' 



Hubble 

type /3 



Disk component 



M./L (Bol-V) (U-V) (B-V) (V-K) 



Whole galaxyC") 
A//L 



Sa 


0.8 


2.42 


-0.78 


0.60 


0.63 


2.70 


0.033 


4.48 


Sb 


0.5 


1.62 


-0.78 


0.44 


0.56 


2.59 


0.030 


3.06 


Sc 


0.1 


1.10 


-0.79 


0.34 


0.51 


2.51 


0.028 


1.86 


Sd 


-0.3 


0.84 


-0.80 


0.28 


0.48 


2.46 


0.027 


1.05 


Im 


-0.8 


0.66 


-0.81 


0.24 


0.46 


2.42 


0.026 


0.66 



From^uzzoin ^2005) template galaxy models 
Including the spheroid component {i.e. bulge + 



halo) 



Stem from the analysis of the theoretical output of our exercise, 
mainly dealing with the expected age-metallicity relation (AMR) 
and with the aged problem of the G-dwarf metallicity distribution. 



3.1 The age-metallicity relation 

By combining eq. ((Sjl and ( 122b with eq. J24b we easily obtain in 
explicit form the expected AMR for different evolutionary scenar- 
ios. With little arithmetic, recalling the £„ and Mo definitions, the 
AMR can eventually be set in terms of the integrated properties of 
the galaxy stellar population: 



0.071 



+ (0.23 + b) (M/L)i 



,0.23 , O.SN 
■ tg + ^ , 



(25) 



providing to express the reference bolometric M/L ratic|f| at 15 Gyr 
in solar unit (as from Table[T), and time tg in GyrQ 

Quite interestingly, note from eq. ( |25b that time evolution of 
Z is almost insensitive to the SFR details, once considering the 
M/L vs. h inverse dependence, as suggested in eq. dlOb . This is 
a not so obvious property of the relation, and derives from the 
nearly constant burning efficiency factor T of Table [T] that con- 
trasts the amount of released energy from mass processing (namely, 
the So/ Mo parameter) with the maximum energy allowed per unit 
mass (i.e. Wt in eq. [8]l. One sees, from Table [T] that F varies by 
less than 20% when moving from a "bursting" star formation event 
to a roughly constant or even increasing SFR with time. This makes 
both Myz and M* to scale with 6 in a similar way. For this reason, 
the same aproximate Z oc t^^'^ dependence has to be expected for 
a wide range of evolutionary scenarios, including the case of SSP 
star bursts. 

The SN contribution is accounted for, in eq. (125b , by the term 
^SN _ + as from eq. ( 122b . where its explicit age de- 

pendence, according to the different evolutionary regimes, is con- 
strained by the A{t, A, 6) function, as in eq. l |23b . As shown in 
Fig. |2l the effect of SN mass processing quickly reaches a steady 
contribution to Z yet within the first Gyr of the galaxy life. For this 
reason, we can therefore safely set the SN enrichment to 



Z^ 



: 0.0066/(1 + R) 



(26) 



° Again, the ratio refers to the stellar mass, that is the total mass ever con- 
verted to stars. 

^ A less direct, but still analytically manageable, form for [Fe/H] vs. t 
could also be obtained, in case, by recalling the basic detinition: [Fe/ H\ = 
\o%Z - log[0.77 - Z{1 + R)] + 1.54, where we assumed (Zq, Yq) = 
(0.02, 0.28) for the Sun, and a primordial Helium abundance Yp = 0.23. 
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Figure 3. Yield chemical abundance for the CSP set of Table 3. Theoretical 
output is from eq. 1241 . assuming a range for the enrichment ratio R = 
/^Y/ l\Z = 3 it 1, according to the label in each panel. Note the nearly 
insensitive response of both y and Z to the SFR details. 



for any practical application in the framework of our discus sionlf] 
After accounting for the eq. ( 110) fit, the AMR of eq. ( I25t reduces 
to 



Z(t) 



0-068 0.23 
1 + 0.18 6 ' 



' 0.0066 



(1 + i?)- 



(27) 



The equation confirms that chemical enrichment of a galaxy disk is 
a potentially fast process. With i? = 3, for instance, the solar metal- 
licity is quickly reached yet within the first few Gyrs of the galaxy 
life, and a current figure about 2 ~ 1.5 Zq after one Hubble-time 
evolution. 

The detailed evolution of yield metallicity along the entire 
Hubble morphological sequence can easily be obtained from the 
corresponding h values, as summarized in Table|2] A more general 
display of the {Z, y) evolution, for CSPs of different SFR power 
index /3, is displayed in Fig.[3]along the assumed range for the en- 
richment ratio R. The relevant case for ii = 3 is also reported in 
detail in Table[3] 



3.1.1 The AMR in the Galaxy 

A summary of the predicted chemical evolution of the Galaxy 
from some of the most recognized models in the recent lit- 
erature is proposed in Fig. |4] In particula r, the m od els by 



iMatteucci & FrancoisI (|l989h. Iwvse & Silk 



1 19941), 



Pardi & Ferrini 



iTimmes. Wooslev, & Weavei* 



Jl994 7Ji 



1995) 



1996) 



iPrantzos & Auberl 
1995^, 




Giovagnoli & Tosi' 
Pilyugin & Edmunds ( 1996), Mihara & Takahara 
as well as by tC hiappini, Matteucci, & Gratton ( 1997) 



IPortinari, Chiosi, & Bressaa (.1998.) , .Boissier & Prantzos. (,199ft) . 
and lAlibes. Labav, & Canal ( l200lh are considered, along with our 
Z law for 6 = 0.5 ± 0.5, from eq. ( |27] |. 

The whole set of theoretical AMRs is compared in the 



* To a closer analysis of Fig.|2] our aproximation tends to slightly overes- 
timate the right figure for more "bursting" SFRs, as in the case of Sa 
galaxies. This has no larger impact, however, on our conclusions, as even 
a fully crude lower estimate such as Z^^ ~ 0.0047/ (1 + R) would lead 
current [Fe/H] predictions for these galaxies to change by much less than 
0.03 dex. See Sec. 4 for an appropriate context to assess this issue. 



Table 3. Yield chemical abundance {Z 
and enrichment ratio R = (Ay/AZ) = 



y) for CSPs with different SFR 

3(t) 



R=3 



Flat or decreasing SFR with time 



tg 


0.8 


0.6 


04 


0.2 


0.0 


[uyrj 




(z . y) 


(z , y) 


(z . y) 


. J'} 


1 


0.019 0.286 


0.017 0.281 


0.016 0.279 


0.016 0.278 


0.016 0.277 


1.5 


0.020 0.290 


0.019 0.286 


0.018 0.283 


0.018 0.282 


0.017 0.281 


2 


0.021 0.294 


0.020 0.289 


0.019 0.287 


0.019 0.286 


0.018 0.285 


3 


0.023 0.300 


0.022 0.295 


0.021 0.292 


0.020 0.291 


0.020 0.289 


4 


0.025 0.305 


0.023 0.299 


0.022 0.296 


0.022 0.295 


0.021 0.293 


5 


0.026 0.308 


0.024 0.302 


0.023 0.299 


0.023 0.298 


0.022 0.296 


6 


0.027 0.311 


0.025 0.305 


0.024 0.302 


0.023 0.300 


0.023 0.299 


8 


0.029 0.317 


0.027 0.310 


0.025 0.306 


0.025 0.305 


0.024 0.303 


10 


0.030 0.321 


0.028 0.314 


0.027 0.310 


0.026 0.309 


0.026 0.307 


12.5 


0.032 0.326 


0.029 0.318 


0.028 0.314 


0.028 0.312 


0.027 0.311 


15 


0.033 0.330 


0.031 0.322 


0.029 0.317 


0.029 0.316 


0.028 0.314 


R=3 




Increasing SFR with time 
















tg 


-0.2 


-0.4 


-0.6 


-0.8 


-1.0 


[Gyr] 


(z . y) 


(Z . y) 


(z . y) 


(z . y) 


(z . y) 


1 


0.016 0.276 


0.015 0.276 


0.015 0.275 


0.015 0.274 


0.015 0.274 


1.5 


0.017 0.281 


0.017 0.280 


0.016 0.279 


0.016 0.278 


0.016 0.278 


2 


0.018 0.284 


0.018 0.283 


0.017 0.282 


0.017 0.281 


0.017 0.281 


3 


0.020 0.288 


0.019 0.287 


0.019 0.286 


0.019 0.286 


0.018 0.285 


4 


0.021 0.292 


0.020 0.291 


0.020 0.290 


0.020 0.289 


0.020 0.289 


5 


0.022 0.295 


0.021 0.294 


0.021 0.293 


0.021 0.292 


0.021 0.291 


6 


0.023 0.298 


0.022 0.296 


0.022 0.295 


0.022 0.295 


0.021 0.294 


8 


0.024 0.302 


0.024 0.301 


0.023 0.300 


0.023 0.299 


0.023 0.298 


10 


0.025 0.306 


0.025 0.304 


0.024 0.303 


0.024 0.302 


0.024 0.301 


12.5 


0.026 0.309 


0.026 0.308 


0.026 0.307 


0.025 0.306 


0.025 0.305 


15 


0.028 0.313 


0.027 0.311 


0.027 0.310 


0.026 0.309 


0.026 0.308 



^ The values of (Z' , ' ) for any other enrichment ratio R' can be derived from the R ■■ 
case, as Z' = [4/(1 + R' )] Znniy' = y + [{R' - 3)/{R' + 1)]Z. 



figure with the empirical AMR in the solar neighborhood, as 
obtained by different stellar samples, probing the young O- 
B stars and the unevolved F-G subgiant population. In par- 
ticular, six important contributions still provide the reference 



framework for this a nalysi s, namely the work bv |Twarog 



1980l) 



199ll) . 



Carlberg et al.l ( 1985h, iMeusinger, Stecklum. & Reimanr 



Edvardsson etal] ( Il993h , iRocha-Pinto et alj ( I2OO0I) , and 



the exhaustive Geneva-Kopenhagen catalog of stars in the so- 
lar neighbourhood (NoMst,rojm et a l. 2004), recently revised by 
iHolmberg, Nordstrom. & Andersenl hoOlh . 

Just a glance to Fig. |4] makes clear the wide spread of the- 
oretical predictions, that span nearly a factor of 3 (namely 0.3- 
0.5 dex in [Fe/H] range) in the expected Galaxy metallicity at a 
given epoch. For example, newly bom stars at f = 10 Gyr yet 
exceed the solar value in the IMatteucci & FrancoisI d 19891) model, 
being as rich as \Fe/H] = +0.11; this is roughly twice the value 
of the IPortinari. Chiosi. & BressanI (|T998) model, which predicts 
on the contrary [Fe/H] = —0.13. Instead of any further specula- 
tive attempt to elect the "right" solution amon g the many envisaged 
scenarios (see, however. iRomano et al.l20()5l for an interesting dis- 
cussion on this issue), we could rather look at the different theoret- 
ical outputs of Fig.|4]as a plain evidence of the inherent uncertainty 
of any "deterministic" approach to Galaxy chemical evolution. 

It may be useful, however, to try a comparison between mod- 
els and observations in terms of the same fitting parameters, by 
aproximating the different AMRs with a general power-law evolu- 
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Figure 4. The solar-n eighbourhood AM R, according to different observin g surveys, as labelled in the plot. I n particular, the works of[Tw aro° (1980') (dots), 
ICarlberg et"al]h985l) (•triangles). .Meusinger. Stecklum. & ReimannHl99lh fsguaresl lEdvardsson et al.Nl993l) (diamonds), and lRocha-Pinto et al.(.200d ) (pen- 
tagons) have been considered. Observations are compared with model predictions from several theoretical codes (the shelf of thin curves; see Fig. |5]for a 
detailed source list). The expected evolution of yield metallicity, according to eq. Ml\ . is also displayed, assuming {Fe/ H] = log Z/Zq, and a birthrate 
b = 0.5 ± 0.5, as labelled on the curve. A current age of 13 Gyr is adopted for the Milky Way. See text for a full discussion. 



tion such as Z oc t'' or, equivalently, [Fe/H] — ("logtg + 1^0 In 
this notation, ( is the slope of the [Fe/H] vs. log t relationship, as 
derived from the fit, while [Fe/H]io = + uj provides the nomi- 
nal metal abundance of newly born stars at i = 10 Gyr. Table|4]and 
Fig. [5] summarize the results of our fit for the whole set of mod- 
els analyzed, comparing with the corresponding parameters from 
eq. (|27ll for i? = 3 and fe = 0.5 ± 0.5. 

As a striking feature, the plot clearly shows the stronger chem- 
ical evolution generally implied by the theoretical codes. The lat- 
ter predict, in fact, a steeper slope for the [Fe/H] vs. logt rela- 
tion, exceeding the observed [Fe/H]io metallicity and predicting 
a larger value of i^. To some extent, this untuned behaviour adds 
further arguments to the delicate role of the gas-infall mechanisms, 
so extensively invoked by modem theory. Actually, while sustain- 
ing star formation, any external gas supply may even ease metal 
enhancement, rather than moder ate it by dilu t ing processed mass 
into a larger amount of fresh gas (lLarsorJl976l : lFirmani & Tutukovl 
Il992h . 



3.1.2 Constraints to the Galaxy SFR 

For its physical definition, yield metallicity provides an upper en- 
velope to the observed AMR, throughout (see, again. Fig. |5j- In 
this regard, it may be useful to further explore this concept, in or- 
der to better clarify the link between implied Z evolution and the 
observed galaxy AMR. 

If we assume that a net fraction f{t) of the stellar mass is 



^ Of course, this aproximation holds as far as the SN term could be ne- 
glected, that is for Z{t) > Z^'^ . 



Table 4. Fitting parameters for Galaxy theoretical and empirical AMR.'"' 





Yield metallicity 




Source 


C 


[Fe/H]io 


This work*''' 


0.23 


+ 0.17 ±0.04 




Models 




Source 


C 


[Fe/H]io 



^atteucc^^raajois Q98^ 
Wvse&Silk(1989> (n=l) 
(n=2) 

Cariei (1994) 
.Paidi&Ferrini(1994) 
Prunlzos& Aubert(1995) 
Timraes, Wooslev. & Weaver (1995) 
Giovagnoli & Tosi (1995) 
Pilvusiii & Edmunds (1996) 
Mihara & Takahara (1996) 
Chiappini, Matleucci, & Gratton ( 1997) 
Porlinari. Chiosi, & B ressan ( 1998) 
Boissier & Pranlzos ( 1999) 
Alibes. Labav. & Canal (2001) 



0.84 
0.79 
0.49 
1.14 
1.04 
0.76 
0.42 
1.06 
0.81 
1.28 
0.71 
0.73 
0.82 
0.73 



+0.11 
-0.07 
-0.04 
+0.08 
+0.00 
+0.00 
+0.02 
+0.09 
-0.06 
-0.00 
-0.00 
-0.13 
+0.02 
-0.03 



Milky Way 13 Gyr old 


Observations 

C 


[Femiio 


Twaroi!(1980) 




0.63 ± 0.04 


-0.01 ± 0.03 


Carlbereetal. (1985) 




0.56 ± 0.08 


+0.02 ± 0.07 


Meusineer, Stecklum. & Reimann 


1991) 


0.54 ± 0.06 


+0.02 ± 0.05 


Edvardsson etal. (1993) 




0.62 ± 0.12 


-0.10 ± 0.11 


Rocha-Pinto etal. (2000) 




0.88 ± 0.04 


-0.05 ± 0.04 


Holmberg, Nordstrom, & Andersen 


(2007) 


0.30 ± 0.01 


-0.12 ± 0.01 


Mean of all surveys 




0.59 ± 0.08 


0.04 ±0.03 



'"■^ Assuming [Fe/H] — log tg + u;, with tg in Gyr. By definition, [Fe/H] iq — -\- u 

is the expected metallicity at i — 10 Gyr. 
'''^ From eq. |25| , for R — 3. [Fe/ H]iq range limits from varying galaxy birthrate b — 
0.5 ±0.5. 
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Figure 5. Representative parameters for the observed (dots) and predicted 
(squares) AMR for the sol ar neighborhood. The theoretic al works by 
iMatteucci & FrancoisI il989l) (labelled as "MF" on the plot). IWvse & Silkl 
(WSl and WS2, respectively for a Schmidt law with n = 1 and 2) , 
Lcaiigi ( 1994) (C), Pardi & Fenini ( 1994) (PF), Prantzos & Aube rJ lll995h 
(PA) iTimmes W ooslev. & Weaver 11995) (TWW), Giovagnoli & Tosil 
jl995|) (GT). [Pilvugin & Edmunds (1996) (PE), Mihara & Ta kaharai 
Il99a) (MT), jO iiappini, Matteucci, & Gratton "Tlgg^ ( CMG) , 
IPortinari. Chiosi. & Bressan H99&) (PCB ), Boissier & PrantzosI h999h 
(BP), and lAlibes. Labav. & Canall <200ll) (ALC) have been acounted 
for, including the expected output from the Z law of eq. )27t (square 
marker labelled as "B"). The Milky Way is assumed to be 13 Gyr old 
for every model. The theoretical data set is also compared with the 
observations, a c cordin g to the surveys of Twarog ( 1980) (dot number "1") , 
ICarlb erg et aP ('1985') f "2"). iMeusing er. Stecklum. & Reimann <199ll) 
("3" ). Edvardsson et ai] h993l) ("4"), .Rocha-Pinto et al. {2Q0Q) ("5"), 
and iHolmberg. Nordstrom. & AndersenI boOVl) ("6"). An AMR in the 
form [Fe/H] = f logtg + uj, is assumed, throughout. The displayed 
quantities are therefore the slope coefficient f , vs. the expected metallicity 
at t = 10 Gyr, namely [Fe/H]io = + ui. Note that models tend, 
on average, to predict a sharper chemical evolution with respect to the 
observations. This leads, in most cases, to a steeper slope and a slightly 
higher value for [Fe/H]io. See the text for a discussion of the implied 
evolutionary scenario. 



As expected, eq. l |29l l shows that, when fresh gas vanishes (i.e. Q — >■ 
0), then Zgas tends to Z. 

According to eq. ( |28t , the total (i.e. fresh + processed) gas 
fraction can eventually be written as 



gtot = G + fs = g{i- f) + f. 



(30) 



Whatever the specific model assumptions about the mass recycling 
mechanisms, it is clear from the equation that / ends up by driving 
C/tot when the fresh-gas reservoir is going to be exhausted. This 
happens, of course, at the very late stages of galaxy evolution or. 
Relying on this definition, eq. ( |29l > can also be conveniently ar- 
ranged in an alternative form such as 



Z 



f 



f 



l-Gt. 



(31) 



As, by definition, / ^ 1 a further interesting property of 
eq. i29\ is that (Zg^s/Z) ^ S. We can take advantage of this sim- 
ple constraint to set two important boundary conditions to galaxy 
SFR. First of all, as 5 cx t^"'', then it must be 



log Zgas - log .Z > (1 - /?) log t + const. 



(32) 



The l.h. side of the relation (or its equivalent form in terms of 
[Fe/H]) simply contrasts the observed AMR in the solar neighbour- 
hood with the expected "yield" metallicity evolution of eq. l l25b . 
The corresponding analytical functions, for the Milky Way case, 
can be found in Table |4] In particular, by comparying the relevant 
slopes, C^, of the [Fe/H] vs. log / relations one eventually concludes 
that ((AMR ^yieid^ ^ (1-^5), which implies, for a current age 
of 13 Gyr, that 

b ^ (1 - 13) < OA ±0.1. (33) 

The latter constraint simply derives from the fact that, evi- 
dently, (1 — g), as well. Therefore, at any epoch, 
the mass fraction of fresh gas in the galaxy could at most be 



e < 1 - 



z 



(34) 



Again, by relying on the fitting functions of Table |4l for a 13 Gyr 
old Milky Way one sees that log(Zgaa/.2) ~ [Fe/H]AMR — 
[Fe/i/]yioid = -0.17 ± 0.05 + log[{l + 7?)/4], so that the cur- 
rent mass fraction of primordial gas in the Galaxy (Qo) must locate 
within the range 



returned to the ISM within time i I by polluting the residual com- 
ponent of fresh primordial gas (Mgas), then the resulting ISM (or, 
shortly, the "gas") metallicity, ^gas, must be 



fZ 



(Aftot/Af*) + (/ 



(28) 



This quantity can be regarded as the representative metallicity of 
newly bom stars. If <S = Af«/Aftot is the mass fraction of "pro- 
cessed" mass and, accordingly, G = 1 ~ S the fresh (i.e. "unpro- 
cessed") gas fraction then, with little arithmetic, eq. ( I28t reduces 
to 



Z 



+ 1 



-1 


'1-5 / 




[ .fs +1 



(29) 



Note that fraction /(t) has to be intended as the net balance of the 
global amount of processed mass retumed to the ISM minus the fraction 
of enriched mass engaged again in later star formation processes. 



OsC^o^ 0.17 (4.9 -_R). 

± 8 



(35) 



By the way, this also implies that R < 5, thus setting an useful 
constraint to the Galaxy enrichment ratio, too. 

Overall, once compared with the representative SFR proper- 
ties of the other spiral galaxies alo ng the Hubble morphological 
sequence teoberts & HavnesI 1 1994I) . these figures are consistent 
with an intermediate Sb type for the Milky W a y (e.g . lvan der KruiJ 
Il986h . In particular, according to iBuzzon i (l2005h . the allowed 
range for the SFR leads to an integrated disk color abou t (B — V) ~ 
0.60±0.05, in agreement with the emipir ical evidence jFlvnn et al.l 
l2006l : lMelchior. Combes. & Gouldl2007l) . 



3.1.3 The returned mass fraction 

The exact knowledge of the way processed stellar mass is retumed 
to the ISM, and the disaggregated contribution of the different ele- 
mental species is clearly the real core of any theoretical effort aim- 
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Figure 6. The expected fraction -FgSP = Mret / Mtot of stellar mass returned to the ISM for a set of SSPs, with changing the IMF power- law index (s), a s 
labelled to the right. The evolution is traced both vs. SSP age and the main-sequence Turn Off stellar mass (mxo). the latter via eq. (3) of lBuzzoii 
The relevant case for a Salpeter IMF (s = 2.35) and two indicative variants for giant- (s = 1.35) and dwarf-dominated (s = 3.35) SSPs are singled out in 
both panels, while the thin dotted curves account for the intermediate cases, by steps of As = 0.05 in the IMF index. 



ing at tracing the physical details of galaxy chemical evolution. In 
our straight approach, this concept is fully condensed into returned 
mass function f{t), which is the ultimate key to quantitatively 
assess eq. \29h The function includes, in principle, a formidable 
wealth of input physics, dealing with chemical yields for individ- 
ual stars of different mass, the quiescent and explosive mass-loss 
properties of stars, the galaxy vs. external environment interplay 
etc. Such a wide range of ingredients, and their related uncertain- 
ties, gives actually reason of the remarkable differences among the 
Milky Way models reviewed in previous section. 

A very plain, but still instructive approach to the problem 
of ISM pollution by returned stella r mass has been investigated 
by iKennicutt. Tamblvn. & Congdonl j 19941 ). following the classi- 
cal scheme of the instantaneous gas recycling extensively prac- 
ticed in the past literature (see, e.g. iLarson. Tinslev. & Caldwell 
Il980l) . Although "old-fashioned" and evidently inadequate for 
any detailed assessment of specific elemental abundances (e.g. 
lArim oto & Yoshii 1986), this method may however still accomo- 
date, although with some important refinements, within our analy- 
sis for setting a safe upper li mit to /. Following the classical treat- 
ment (see, e.g. lTinslevlll980h . we can consider a SSP of total mass 
Mtot, and an IMF stellar mass range between 0.1 and 120 Mq. 
With these constraints, the fraction of mass returned to the ISM 
within the lifetime of Turn Off stars, of mass mro, quantifies in 

M,et {2~s) fi20 
^^^"^ = ir~ " [^2-^1120 / A™(m)m dm. (36) 

Mtot [m Joi Jm{TO) 

This estimate requires to know the Am(m-) function, that is 
the mass amount lost by stars of initial mass m, during their 
Post-MS evolution, so that Am , = r n — mfln. Following the 
IKennicutt. Tamblvn. & Congdonl (1993) arguments, we could rely 
on fully empirical facts, and take the observed Am rela tion, as from 
the Ga lactic open clusters. For example, according to IWeidemannI 
a nice fit can be obtained for stars within the range 
0.5 Mq < m < 8 Mq, such as 

A„ = 0.916 m- 0.44 (37) 
±4 ±2 



For higher masses, if we assume that a Chandrasekhar core of 
1.4 Mq is left all the way by the SN events, then 

A„ = (m-1.4). (38) 

for m > 8 Mq . With these figures, the results of our calculations 
are shown in Fig.|6l exploring the trend of Fssp with varying SSP 
age or, equivalently the stellar mass marking the MS Turn Off point 
of the population. Note that Fssp smoothly increases with SSP age, 
and for a Salpeter IMF it reaches a "ceiling" value about Fssp = 
0.32 after one Hubble timeEl 

On the other hand, as far as a CSP with a more entangled star- 
formation history is concerned, we have to expect the total fraction 
of returned mass, Fcsp{t), to be the result of a convolution with 
the SFR, such as 

FcsAt)^ [ Fssp(r) ® SFR(r)dr. (39) 
Jo 

This leads, in general, to a lower value for -Fcsp(t), compared 
to the corresponding SSP case. In addition, for a proper assess- 
ment of the gas budget in the galaxy one has also to consider that 
part of the returned mass is in fact to be recycled again into fresh 
stars, its exact amount being indeed a crucial distinctive property 
of any theoretical scheme aimed at tracing galaxy chemical evo- 
lution. Our conclusion is therefore that, in any case, it must be 
/(i) =SFcsp(t) <Fssp(t) <0.3. 

3.2 Mixing processes and metallicity spread 

One cominon feature of all these models, however, is that f{t), and 
the physical conditions that led to establish the galaxy AMR are 
univocally defined at any epoch. As a consequence, by definition, 

One has also to be aware, however, of the extreme sensitivity of this up- 
per Hmit to the IMF power index, as evident from the curve grid in Fig.|6] 
A relevant case in t his regard is that of a "dief Salpeter IMF, as paramet er- 
ized for instance bv lKennicutil il983h or lKroupa. Tout. & Gilmoi3 h993l) to 
better account for the observed flattening in the dwarf-star number counts at 
sub-solar mass. By means of eq. j36t one could verify that this IMF closely 
traces the standard Salpeter case being however capped at a slightly higher 
value of Fssp ~ 0.46. 
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0.2 0.4 0.6 0.8 1 

Mass fraction of fresh gas 

Figure 7. The output of eq. )29t is displayed vs. the mass fraction of fresh 
primordial gas (Q). The two shelves of curves assume the fraction / of 
returned stellar mass as a free parameter. The lower shelf (dashed curves) 
is for / = 0.01 0.09 by steps of A/ = 0.02, while the upper shelf of 
curves (soHd lines) cover the / range from 0.1 to 0.3, by steps of A/ = 0.1. 
The shaded area within —0.3 ^ log(Zgas/.Z) ^ single out the allowed 
f-Q combinations that may account for the observed metallicity spread of 
stars in the solar neighborhood. 



no spread in [Fe/H] is envisaged by theory among coeval stars, 
throughout. On the other hand, observations of unevolved nearby 
stars clearly point to some degree of metallicity dispersion, perhaps 
slightly amplified by observational unc ertainties and the adopte d 
age/metallicity classification procedure (Rocha-Pinto et alJ IlOOOh . 
For example, by relying on the Edvardsson et al. ( 1993) survey, we 
derive (T[i^e/ff] = 0.23 dex from the distribution of the individual 
stars around the mean AM R locus. This i s about twice the typi- 
cal value originally found bv lTwaro^ d 1980l) . an d can be contrasted 
with cj \Fe/H] = 0.16 dex, as obtained for the lRocha-Pinto et al.l 
(I2OOOI) stars, or with a\Fe/H] = 0.18 dex, displayed by the 
iHolmberg. Nordstrom. & Anderse 

3 (l2007b extended star sample. 
Whether this spread is really the gener al case for spiral galaxies is 
a still open and debated qu estion (e.g. iKoeppen & Arimotolll990l : 
IWvsell2006l : |Prantzosll2008h . 

Again, the envisaged Z evolution can help, we believe, to set 
some interesting constraints to this problem, as well. Following 
eq. ([29}, Fig. [7] summarizes the expected response of the Z^^j Z 
ratio to just the two reference parameters that modulate galaxy 
chemical evolution in our framework, namely (J, and /. When re- 
garded into a broader context, our plot displays expected relationa- 
ship of the global mass fraction of fresh gas in the galaxy. Alterna- 
tively, this "macro" view can be complemented by a "micro" view, 
being the Q parameter intended as the local gas density of star- 
forming regions, such as t/ = Pgas/ptot, and featuring therefore 
the specific mixing conditions that embed star formation. One sees 
from the figure that, in principle, a wide range (a few dex!) may 
be allowed to ^gas, at every epochs, as coeval stars could naturally 
spread in the plot along an / = const locus, in consequence of the 
micro-environment conditions proper to the different star-forming 
regions i nside a galaxy. 



The lidvardsson et al.l (Il993[). iRocha-Pinto et"ai] ( I2OO0I) , and 
iHohnberg. Nordstrom. & Anderse 3|(l2007') star samples offer once 
more an important feedback in this regard, as shown in Fig. |8] For 
the first sample, 161 stars are nominally younger than 13 Gyr, the 
assumed age of the Milky Way; the overall statistics of their indi- 
vidual A [Fe/Jf] = [Fe/ii"]yij.id-[Fe/Ji"], distribution indicates 
that 88% of the stars are comprised within A[i^e/ff] ^ 0.50 dex. 



while a fraction of 65% lie within A[i^e/ff] 0.30 dex. With the 
same p rocedure, for the 412 relevant stars in the iRocha-Pinto et al.l 
I2OOO) sample we obtain a fraction of 97% and 84%, respectively, 
within the same A \Fe / H] limits. Similar figures (namely 97% and 
83%, respectively) apply also for the 13141 eligible stars in the 
Geneva-Kopenhagen calatog. Allover, this means that over 4/5 of 
all newly-born stars in the solar vicinity never departed by more 
than 0.3 dex from Z, at any epoch. 

Facing a so tight A[Fe/H] = log(2/Zgas) star distribu- 
tionH one may actually wonder why stars in the Galaxy appear 
so closely "timed" with the yield metallicity, everytime. 

In spite of any more or less sophisticated approach to this cru- 
cial issue, just a glance to Fig. [7] makes evident that the observed 
range for the Z^^^jZ ratio univocally implies (at least in the so- 
lar vicinity) that star formation must have proceeded within ISM 
conditions characterized, all the way, by a scanty gas abundance 
(G ^ 0.5, see again Fig.^Tj). 

Overall, this feature seems to add further strength to well 
recognized non-standard scenarios, where the plain depend ence 
of SFR on gas density, like in a classical ISchmidtl jl959l) law, 
might more effectively be replaced by some positive feed-back 
to other triggering mechanisms, either related to disk dynami- 
cal properties, inefficient gas mixing or molecular cloud bunchin 



jLarson & Tinslevll978l:lwys e & Sil kll"98^;lDoDita & Rvdeill99- 
Oey. 200Q; 'Soissier et alj"| So3 : Scal o & Elmegreenll2004l: iNittlej 
I2OO5I : Ivan Zee & Havn es'%06 ) . 

A further interesting piece of information can be added to the 
emerging scenario when probing the A[Fe/H] distribution at dif- 
ferent Gal actocentric distance. This can b e done, for instance, by 
relying on Carraro. Ng. & Portinaril(ll998b . who collected homoge- 
nous data for a set of 37 open clusters. The evident advantage, in 
this case, is that age and metallcity can be derived in a more direct 
way, based on the observed c-m diagram of the whole stella r mem- 
bers of each clusterF^ The lCarraro. Ng. & Portinaril(ll998h results 
are reported in the lower panel of Fig. |8] marking the cluster dis- 
tance from Galaxy center with a different dot size (i.e. the closer 
the cluster, the bigger its marker). 

The local properties of disk metallicity are consistently traced 
also by the open cluster population, whith a clear additional evi- 
dence, however, for less pronounced chemical enhancement pro- 
cesses in the most peripheral stellar systems, for which a systemat- 
ically larger departure from Z is displayed in Fig. [8] 



3.2.1 The "G-dwarf" problem 

In its essence, the tuned distribution of the A[Fe/H] parame- 
ter, as in previous discussion, has much to do with the classical 
"G-dwa rf" problem (van den Bergh 1962; Pagel & Patchett 1975j; 
lTwarodrr98 0), that is an evident lack, nowadays, in the [Fe/H] 
(or log Z) domain of the "desired" population of extremely metal- 
deficient unevolved stars, reminiscent of the earliest evolutionary 
phases of the disk formation in the Galaxy. 



Strictly speaking, for its delayed supply to the ISM, Iron may not be a 
fair proxy of the whole metallicity among very metal-poor (i.e. "pristine") 
stars in the Galaxy (see the following Sec. 3.2.2 for a most explicit assess- 
ment of this issue). This assumption canies therefore some uncertanties for 
stars within t < 1 Gyr in the plots of Fig. [8] 

One has to remark, however, the evident bias of the 
ICarraro. Ng. & Portinaril ( Il998l) sample against very old open clusters 
likely due, according to authors' discussion, to the disrupt ive mechanisms 
along the past dynamical history of the Galaxy jFrieil99a) . 
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Figure 8. The observed metallicity distribution of stars in the solar 
neighbourhood. The displayed quantity is A[Fe/H] = [Fe/H], — 
[Fe/ H]yigi^ = log(Zp;as/2)- Different stellar samples have been consid- 
ered froni_the_woA_ofjEdvardsson et al. ( 1 993 ), Rocha-Pinto et al. 1 20(3). 
and lHolmberg. Nordstrom, & And ersen (2007), from top to bottom, as la- 
belled in each panel. The low er panel dis pla ys the sam e trend for a sample 
of open stellar clusters, after ICarraro, Ng, & Portinaril (Il998l) . The differ- 
ent Galacticentric distance for these systems is indicatively marked by the 
dot size (i.e. the smallest the dot, the largest the distance from the Galaxy 
center). For each distribution, the relative fraction of stars within a -0.3 and 
-0.5 dex A[Fe/H] residual is reported in each panel. A cumulative his- 
togram of the data is also displayed, along the right axis, in ai'bitrary linear 
units. 



According to our SFR parameterization, for 6 < 0.4, an ex- 
pected fraction of ~ (50 10^15 10^)'' > 10% of G-dwaif stars 
should have produced in the disk of our galaxy within the lapse 
(t < 50 Myr) of the first type-Il SN burst. According to Fig. |2), 
for these stars a metallicity as poor as [Fe/H] ~ —1.10 should be 
expected. The lack, nowadays, of such a relevant population of vir- 
tually "zero-metal" stars is actually the central issue that led theo- 
rists in the late 80's to invoke a composite formation for the Galaxy 
disk, through a "buffering" dynamical mechanism that led first to 
the forma tion of an outer "thin" disk, followed by a n inner"thick" 
structure ( lOilmore & Wvselll986l:lFerrini et al.lll992l) . 

Translated into a linear metallicity scale, the undersized pop- 
ulation of stars with [Fe/H] <C — 1 can be matched by a roughly 
uniform star density along the Z range, so that dN,/dZ ~ 
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Figure 9. The [Fe/H] distribution of G-dwarf s tellar samples in the solar 
neighborhood (shaded histograms), according to iRocha-Pinto et al] j2000l) 
(upper panel), and Nordstrom et al. (2004) (lower panel), is compared with 
the expected yield-metallicity distribution accounting for the different star- 
formation histories as for the disk stellar population of Buz zoni ( 2005.) tem- 
plate galaxies along the Hubble sequence (solid curves, as labelled on the 
plot). A roughly constant stellar distribution in the Z domain leads to an 
implied birthrate of the order of fe ~ 0.6, as pertinent to Sbc galaxies (see 
Table 2 for a reference). 



co nst (see, in this regard, Pagell 0*9891 a nd especially the work 
oflBeers. Preston. & Shectmanl ll985. 1992: JSommer-Larsenlll99 ll : 
iRvan. Norris, & Beerslll996l and lCavrel,1996l) ! 

The stellar metallicity distribution can easily be related with 
the other distinctive properties of galaxy evolution in our frame- 
work, as 



dN4Z) dN,{t) 



dZ 



dt 



dZ 
'dt' 



(40) 



The first term in the rh. side of the equation is, of course, propor- 
tional to the SFR (oc t^^), while the second term is simply the 
time derivative of the observed AMR (oc t'""^). By replacing the 
relevant relations, we have 



dZ 



oc Z 



(i-C-/3)/C 



dN, 



cc 10 



[Fe/H-]{l-/3)/C 



(41) 



(42) 



d[Fe/H] 

if one better likes to put the relation in terms of [Fe/H]. For a 
roughly constant dN^/dZ distribution, it must be (1 — — /?) ~ 0, 
or ~ 6. According to Table|4}, for a current age of 13 Gyr, this 
implies that G-dwarf distribution may be reproduced by a Galaxy 
birthrate of 



&~ 0.6 ±0.1. 



(43) 
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A graphical sketch of our conclusions is proposed in Fig. |9] 
by comparing the G-d warf stella r sample of the Galaxy disk by 
iRocha-Pinto & Macie] \\99&) and Nordstrom et al. (2004) with the 
expected metallicity distribution as obtained by applying eq. l l42t 
for a range of SFR representative of the spiral-galaxy sequence in 
theH uzzon il ( l2005h template galaxy models. Note that Milky Way 
observations closely match the standard case of the Sb/Sc morpho- 
logical types (see Table|2}. 

To a better analysis, however, we have to remark the somewhat 
marginal value of 6 when compared to our previous conclusions, 
as the tight relationship between Galaxy AMR and Z trend rather 
called for a lower figure (i.e. b < 0.4). This apparent dichotomy 
actually summarizes the full essence of the G-dwarf problem em- 
phasizing, from one hand, the central role of the early stellar com- 
ponent in constraining the Galaxy AMR, but requiring, on the other 
hand, for these stars to be bom "elsewhere" (i.e. in the halo or in 
the "thin" disk). 



3.2.2 The "a" -element versus Iron-peak enrichment 

A further implication of the simple energetic arguments outlined 
in Sec. 2 is that metal enrichment should have proceeded selec- 
tively in the very early stages of Galaxy evolution, according to 
the composite inteiplay between "quiescent" luminosity evolution 
of stars and violent action of SNe. In fact, bright stars are to be 
regarded as the elective suppliers of lighter pre-Carbon elements, 
while SNell further add to the synthesis of the a-element chain, es- 
pecially providing Oxygen and heavier metals (Woosley & Weaver 
1 199^ . For their different physical conditions, SNela are meant, on 
the contr ary, to be the prevailing donors of Fe a nd other Iron -peak 
elements jNomoto. Thielemann. & Yokoill984h to the galaxy ISM. 

Given the delayed appearence of SNia, one has to expect 
a elements (and Oxygen, in particular) to cumulate earlier than 
Iron-peak elements into ISM of the Milky Way. This feature di- 
rectly calls, therefore, for some 0-Fe decoupling among the low- 
metallicity stellar component of our galaxy. A positive [O/Fe] 
relative abundance is actually a widely recognized property when 
comparing halo and disk stars, and the effect is a supplementary 
piece of evidence dealing with the G-dwarf problem previously dis- 
cussed. The study of the [0/Fe\ vs. [Fe/H] relationship in the 
local framework prove therefore to set important constraints to the 
different feeding channels that eventually led to Z enrichment in 
the Galaxy. 

We can further explore this issue by relying on our model to 
compare the expected [O/Fe] evol ution with the observed trend 
for the Galaxy stars. The work of IPagel & Tautvaisiend d 19951) 
and IScannapiec o & Bildsten ( 2005) can be taken as useful ref- 
erences to set the relevant quantities for Fe and O yields in 
Type la and II SNe. Accordingly, we assume a typical abun- 
dance of [0/Fe\ii = -(-0.50 dex for the released mass by 
SNell and [0/Fe]i = —1.51 dex for the SNia case, being 
log(O/-F'e)0 = -1-1. 33 dex the reference figur e (in number abun- 
dance) for the Sun dGrevesse & Sauv dl ll998h . In addition, with 



IScannapieco & BildstenI 1 2005h . we take a ratio r — {m'p^ 



(0.062 : 0.74) for the released Iron yields according 
to the SN typeEl 



^■^ As Iron is the main output of any SNia explosion, its yield 
must obey an inherent ai'gument of energetic self-consistency being 
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Figure 10. The \0 /Fe] distrib u tion o f stars in the Milky Way is displayed 
acording to the Jonsell et al. (2005) (diamonds) and E dvardsson et alj 
1 1993) (dots) data samples, respectively for the halo and disk stellar popu- 
lations. Our model output, according to eq. | |461 is supeiposed (solid lines), 
for a representative birthrate 6 = 0.4, and for a range of SNia delay time A 
from ^ 50 Myr to ~ 1 Gyr, as discussed in Sec. 2.2. 



Equations ( I13t . ( I16t and i20\ can still accomodate for our cal- 
culations providing to replace the {E51 10^^) / {Wsn Mq) energy- 
release term with the relevant chemical yields. Recalling our SFR 
parameterization, and the notation of Sec. 2.2, the total mass re- 
leased by Type IT and la SNe for the element X can be expressed 



m{X) = m{X)ii + m{X)i B{t, A, fe), 



(44) 



where B{t, A, fe) is a function of age (t) and stellar birthrate (fe), 
and accounts for the SNia delay time (A). Like for eq. ( I23t . pre- 
viously discussed, an analytical form can be devised also for this 
function, providing to express A in years, such as 



B(t,A,fe) 







0.35(1 



for r < 1 



(45) 



l + b ' 

for r ^ 1 



being, again, r — t/A. 

For the [O / Fe] relative abundance we have therefore 



o_ 

Fe 



loa 



. lQlO/Fe]II _,_ ^ -^Q[0/Fe]/ 

r + B 



(46) 



ures of Sec. 2.2.2, we actually have that m'/^ ~ 0.75 ± 0.25 M 



0. 



Note, from eq. (|46}, that B = Ofovt < A, which implies that ISM 
[O I Fe\ abundance simply coincides with the SNiI relative yields at 
the very early evolutionary stages of the Galaxy. Iron enrichement 
is obviously favored by increasing B in the equation. This naturally 
happens when getting the Galaxy older (that is by increasing r) 
and/or by decreasing the SNia delay time A. 

Our results are sketched in Fig. [TO] compari ng with the obser- 
vation of di sk and halo s tars, respectively from lEdvardsson et al.l 
(1993) and Ijonsell eTal] (j2005). In the fi gure, a safe lower and 
upper envelope to the observed [O/Fe] stellar distribution can 
be placed for A ~ 50 Myr and 1 Gyr, respectively. This 
range also consistently comprises the claimed "glitch" in the 
[O/Fe] vs. \FelH\ trend, ideally located about [Fe/ff] ^ -1 



in close agreement with the stan dard model predictions (see, e.g. 
iNomoto, Thielemann. & Yokoll9i4l their model W7). 
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jClegg. Tomkin. & Lamber3ll98lh . Our output therefore confirms 
that most of the SNia impact on the Galaxy metal enrichment 
should have occurred very early within the first Gyr of life, al- 
though just a glance to the data of Fig.[TO]indicates that some 10* yr 
may have required for the SNia process to fully deploy within the 
Galaxy stellar environment. This conclus ion is supported, i ndeed , 
by the empirical arguments of Maoz, S haron,& Gal-Yariil (2010l) 
and th e theoretical ones by IGreggiol (12005 ) and Kobavashi et alj 
( Il998h . 



4 CHEMO-PHOTOMETRIC PROPERTIES ALONG THE 
HUBBLE MORPHOLOGICAL SEQUENCE 

Our quite special lookout inside the Galaxy makes the study of our 
system certainly favored but, at the same time, also subtly biased, 
compared to the analysis of external galaxies. From one hand, in 
fact, we can directly adress the problem of "Galaxy metallicity" 
in terms of the aggregated information from individual stars in the 
solar neighborhood; on the other hand, the criterion to define the 
star sample (either by fixed space volume, or apparent magnitude 
limit, or even selected spectral type etc.) may play a crucial role 
in driving our conclusions, as the "representative" metallicity value 
would eventually depend on the way we account for the contribu- 
tion of individual stars. This is a subtantial difference with respect 
to the study of distant galaxies, for which we usually cannot re- 
solve individual stars at all. In the latter case, one should forcedly 
rely on the integrated spectrophotometric properties of the system, 
as a whole, to derive somewhat "effective" distinctive parameters 
to constrain galaxy properties. 



4.1 Matching the Arimoto & Jablonka (1991) theoretical 
framework 

From the theoretical side, our general view of disk chemical prop- 
erties along the Hubble morphological sequence may usefully com- 
plement the analysis of AJ91, certainly one of the few explicit at- 
tempts in the recent literature to self-consistently tackle the prob- 
lem of chemo-photometric evolution of late-type galaxie s (see, in 
addition, also Arimoto & Yoshii 1986, 1987: Koeppen & Arimotd 
1 19901 to complete the theoretical picture). A first interesting com- 
parison, in this regard, is proposed in Fig. Illlwhere the gas metal- 
licity, log(Zgas/^0), for the 15 Gyr AJ91 disk models along the 
Sa — > Sd sequence is displayed together with the relevant values of 
our yield metallicity, \og{Z / Zq). 

For its different IMF limits (stars are produced between 0.05 
and 60 Mq) the AJ91 stellar populations are slightly "darker" (that 
is with a larger M / L ratio) for their larger fraction of dwarf stars, 
and with a less prominent contribution (~ —20%) of SNe to galaxy 
metal enrichment. These differences with respect to our theoretical 
scheme become more important when moving along the Sa — >■ Sd 
sequence, as stars of increasingly higher mass prevail as contribu- 
tors to galaxy luminosity. 

In spite of the different input physics, one has to remark how- 
ever a notable agreement in the [Fe/H] predictions, as shown in 
Fig. [TT] The case of Sa models is illustrative in this sense, as 
the vanishing residual gas (only 4% of the total disk mass in the 
AJ91 model) makes Zgas to closely approach Z, as expected. Once 
accounting for the increasing fraction of residual gas along the 
morphological galaxy sequence, the prediction of eq. ( 13 U tightly 
matches the AJ91 output. Facing the relatively flat trend of Z 
vs. galaxy type, it is evident from the figure that the increasingly 



0.5 




isi 



O 



-0.5 




2 4 6 8 10 
Morphological Type [T] 



Figure 1 1. The gas metallicity of 15 Gyr disk stellar populations accord- 
ing to the I Arimoto & Jablonkal Jl99lh (AJ91) models (solid line with dia- 
monds) is displayed together with the expected predictions from eq. 1311 
(shaded area) for a range of returned mass fraction / between 0.07 and 0.3, 
as labelled on the plot. For tighter self-consistency in our comparison, the 
adopted gas fraction CJtot for our model sequence matches the con'espond- 
ing figures of AJ9 1 models, while yield metallicity derives from eq. 1271 for 
the relevant birthrate values of the Buzzoni (2005) template models along 
the Sa— >Im morphological sequence (solid line with dots). 



poorer gas metallicity along the Sa Im Hubble sequence is in 
fact mainly the result of the diluting process, that shares enriched 
stellar mass with a larger fraction of residual gas. 

Interestingly enough, in the quite delicate balancing mecha- 
nism that governs the [Fe/H] trend vs. Hubble type, the exact de- 
tails of stellar mass return processes do play a somewhat marginal 
role. This is evident, for instance, by comparing in Fig. II lithe effect 
of changing the return mass fraction / from its allowed maximum 
(/ — 0.3) down to a factor of four lower value (/ = 0.07). Rather, 
the amount of fresh residual gas, that survived past star formation, 
seems a much more central figure to modulate the enrichment pro- 
cess, and this directly calls for the distinctive birthrate b as the ulti- 
mate parameter to constrain disk metallicity. We will further return 
in more detail on this important point in a moment (see Sec. 4.3). 



4.2 "Chemical" vs. "photometric" metallicity estimates 

Due to the coexisting contribution of gas and stars inside a spi- 
ral galaxy, any fair estimate of "galaxy metallicity" i s a far more 
entangled problem. Quoting [Roberts & HavnesI ( 119941) : "intercom- 
parison between abundances based on emission lines from Hll re- 
gions and on absorption lines and colors from a stellar population 
is still uncertain... Atomic processes guide emission-line analysis, 
while the absorption lines must be evaluated in terms of both the 
composite stellar population and the metallicity". 

From the observational point of view, the inherent difficulty in 
disentangling the genuine metallicity effects on the color properties 
of late-type galaxies has greatly favored the alternative approach 
relying on the study o f Hll regions as effective proxies of disk 
metallicity (see, e.g . IPeimbert, Ravo. & Torres-Peimberd Il978l: 
iZaritskv, Kennicutt, &Huchrall 19941 : Peimbert. Carigi. & Peimbert 
[20011 : iGametH l20o3) ~ Broad-band colors and narrow-band 
spectrophotometric indices have been imposing, on the 
contrary, as election too ls to probe chemical abundance 



in the 



as 

galaxy 



19961 



_ bulge s djablonka. Martin. & Arimoto! 

Henrv & Worth evI Il999l: llablonka. Go rgas. & Goudfrooiil 
Gorgas, lablonka, & Goudfrooiil 20071) , while they mainly 



20071: 



Chemo-photometric evolution of spiral galaxies 15 



constrai n galaxy SFR in the disks (among ot h ers, see i n par- 
ticular iKennicutt, Tamblvn. & CongdorJ 1 1994 lOavazzil 1 19931 : 
iBoselli et alJ l200l[) . In this framework, one remarkable effort 
toward a photometric approach to the metall ic ity pro blem is the 
work of iPerez. Sanchez-Blazquez. & Zurital ( |2009|) . providing 
Lick-index measurements of stellar populations in barred spirals. 

An instructive summary of observational galaxy diagnostics 
is proposed in Fig. [121 where the performance of both "chemical" 
(Hll regions) and "photometric" (Lick indices from the aggregated 
stellar population) methods is compared for a set of spirals along 
the full Hubble sequence. To the theoretical pattern of Fig. [TT] re- 
poted as a guide line, we superposed the metal l icity e stimates for the 
galaxy set of Zaritskv, Kennicutt, & Huchra h994') (square mark- 
ers), and for the Perez. Sanchez-Blazquez. & Zurita (.2009) SO/a 
sample (triangles). As usual, the ori ginal gas metallicity in the 
IZaritskv. Kennicutt. & Huchral ( 1 19941) galaxies is given in terms of 
[O I H\ abundance, so that two variants of the figure can be devised, 
depending whether we plainly assume Oxygen to trace global 
metallicity (so that [0/-ff] = [Fe/_ff], as in lower panel of the fig- 
ure) or, perhaps more realistically, whether decoupled Q-element 
enrichment should be considered to rescale chemical abundance 
(upper panel), as we have been discus sing in previous section. In 
the latter case, from different sources ( Clegg. Tomkin. & LambertI 



ll98lHBesseIl, Sutheriand, & Ruanlll99ll : lEdvardsson et al.lll99'3l\ 

weset[0/H] =0.6[Fe///]. 

Compared with the theoretical predictions, it is evident from 
Fig. [12] that any empirical effort to constrain metallicity in spiral 
galaxies still suffers from very high uncertainty. This is especially 
eviden t for the Hll measurements of iZaritskv. Kennicutt. & Huchral 
( Il994l) . where point scatter is actually comparable to the claimed 
internal error of individual measur ements. Ap parently, t his does 
not se em to be the case for the , Perez. Sanchez-B lazquez . & Zurital 
( l2009h Lick data, although the important point-to-point spread casts 
evident doubts on the nominal accuracy claimed for the individ- 
ual [Fe/H] estimates. Definitely, we believe that a a[Fe/H] ~ 
±0.5 dex might be regarded as a realistic figure for current empir- 
ical (either "chemical" or "photometric") estimates of metal abun- 
dance in late-type galaxy disks. 

Once considering in more detail the overlapping range of mor- 
phological types between the two galaxy samples in Fig. [T2| it is 
interesting to note that "photometric" metallicity estimates appear 
to be systematically lower than "chemical" estimates. According to 
our previous arguments, this can be seen as an eloquent example 
of "biased" sampling criteria; in fact, by probing gas metallicity 
we retrieve the current status of galaxy chemical enrichment, while 
through spectrophotometric Lick indices we are probing the galaxy 
stellar population "averaging" along the entire range of age of lu- 
minous stars. 



4.3 Mass, morphology & birthrate (which drives galaxy 
metallicity?) 

Although morphological classification is the most immediate tool 
for a "quick-look" characterization of galaxy properties, it has been 
long questioned whether morphology univocally seats on intrin- 
sic physical properties of the galaxy system as a whole. This is 
especially true for late-type galaxies, where the prominent spiral 
pattern actually involves a relatively small fraction of the galax- 
ian total mass. As extensively discussed in the literature along the 
latest decades jWhitfordll 19621 : lRoberi:s & Havneslll994 iKennicuttI 
Il998l) . the combined analysis of the galaxy leading properties led 
to a number of established correlations between "primary" physical 
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Figur e 12. \Fe/H} metallicity fo r t he t wo galaxy samples 
of j P erez. Sanchez-Blazqu ez. & Zurital j2009l) (triangles) and 
IZaritskv. Kennicutt . & Huchra (squares), compared with the 

yield and gas metallicity figures (solid line with dots and shaded area) as 
from Fig. [TT] While the Perez. Sanchez-Blazauez. & Zurita (2009) esti- 
mates derive from Lick indices tracing the galaxy stellar population, in case 
of Zaritsky. Kennicutt. & Huchra ( 1994) the [Fe/H\ derives from the Oxy- 
gen abu ndance of galaxy Hll regions, by adopting log(O/-ff)0 = 8.83 for 
the Sun iOrevesse & Sauvallll998) . Two cases are envisaged, in this regard, 
assuming that Oxygen straightly traces Iron (and the full metallicity), as in 
lower panel or, more realistically, that a scale conversion does exist such as 
[O/H] = .6 \Fe/ H], as suggested by the observation of disk stars in the 
Milky Way I'Clegg. Tomkin. & Lamberlll98ll : lBe"ssell. Sutherland. & Ruaij 
1 1 99 ll : lEdvardsson et al. 1993). A little random scatter has been added to 
the morphological T class of the points in each panel for better reading. 



param eters, the most popular one being by far the ]Tu11v & Fisheil 
(]l977l. TF) relation. 

On the other hand, the obvious reference role of galaxy lumi- 
nosity, as in the TF relation, can easily be replaced by galaxy stellar 
mass, especially once focussing on galaxy photometric pr operties 
in the mid infrared range, as convincingly demonstrated bv ^Gavazzil 
( Il993h . All these different pieces of evidence led eventually to iden- 
tify galaxy (stellar) mass as (one of) the leading parameters that 
may in timately mark the life and fate of the spiral systems as a 
whole jTullv. Mould, & Aaronsoiilll982l) . 

In this framework, the recent ultraviolet study bv lMarino et al.l 

( 1201 oh may add further interesting arguments to this composite sce- 
nario. These authors collected Galex observations of three spiral- 
rich a ggregates, from the Lyon Groups of Galaxies (LGG; Garci^ 
ll993l) . For their galaxy population, these groups can be regarded 
as close analogs of the Local Group, displaying a quite "clean" 
sample of relatively unperturbed objects seen at comparable dis- 
tances and covering the whole Hubble morphological sequence. 
While Galex data allow us to explicitely assess current SFR in 
these galaxies, the match with long-wavelength magnitudes pro- 
vide, on the contrary, a more direct hint of the amount of stellar 
mass, through the appropriate M/L conversion ratio. 

As a first step of our analysis, galaxy internal extinction has 
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Table 5. Inferred physical properties for the iMarino et al.l fcOlOl) LGG 

galaxy sample 



Galaxy ID 


Morph 

T 

lype 


T 


6 


Mb log Al"*^i 
r r.r. 93 


log SFRo 


b 


NGC 1249 


SBc 


6.0 


-1.02 


-19.57 


10.53 


0.453 


1.07 


NGC 1311 


SBm 


8.8 


-1.07 


-18.30 


9.69 


0.009 


2.67 


IC 1933 


Sc 


6.1 


-1.16 


-18.66 


10.16 


0.226 


1.52 


IC 1954 


SBb 


3.2 


-0.96 


-19.72 


10.93 


0.423 


0.40 


IC 1959 


SBd 


8.4 


-1.29 


-18.21 


9.70 


-0.065 


2.19 








1 nn 1 77 








NGC 1744 


SBcd 


6.7 


-1.18 


-20.24 


10.72 


0.626 


1.04 


NGC 1792 


SBbc 


4.0 


-0.24 


-22.42 


11.91 


1.376 


0.38 


NGC 1800 


Sd 


8.2 


-1.22 


-18.78 


9.95 


0.090 


1.77 


NGC 1808 


SABa 


1.2 


-0.36 


-22.21 


12.16 


1.035 


0.10 


NGC 1827 


SABc 


5.9 


-0.22 


-19.96 


10.70 


0.886 


1.97 


ESO305-009 


SBd 


8.0 


-1.31 


-18.39 


9.82 


-0.022 


1.85 


ESO305-017 


IB 


9.9 


-1.29 


-17.14 


9.10 


-0.504 


3.25 


ESO362-011 


Sbc 


4.2 


-2.85 


-17.74 


10.01 


-1.000 


0.13 


ESO362-019 


SBm 


8.9 


-1.31 


-17.72 


9.45 


-0.314 


2.24 
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NGC 3370 


Sc 


5.1 


-1.24 


-19.76 


10.72 


0.207 


0.40 


NGC 3443 


Scd 


6.6 


-1.18 


-17.26 


9.54 


-0.208 


2.32 


NGC 3447 


Sm 


8.8 


-1.40 


-17.31 


9.29 


-0.014 


6.33 


NGC 3447A 


IB 


9.9 


-1.76 


-15.15 


8.30 


-0.768 


11.0 


NGC 3454 


SBc 


5.5 


-0.62 


-19.08 


10.39 


-0.209 


0.32 


NGC 3455 


SABb 


3.1 


-0.73 


-18.31 


10.37 


0.434 


1.48 


NGC 3501 


Sc 


5.9 


-0.78 


-18.99 


10.31 


-0.142 


0.45 


NGC 3507 


SBb 


3.1 


-1.07 


-20.10 


11.09 


0.422 


0.28 


UGC 6022 


I 


9.9 


-1.31 


-14.87 


8.19 


-1.127 


6.27 


UGC 6035 


IB 


9.9 


-1.27 


-16.48 


8.83 


-0.695 


3.82 


UGC 6083 


Sbc 


4.1 


-0.82 


-17.39 


9.88 


-0.722 


0.32 


UGC 61 12 


Scd 


7.4 


-1.04 


-17.49 


9.54 


-0.248 


2.12 


UGC 6171 


IB 


9.9 


-1.60 


-16.36 


8.78 


-0.847 


3.02 



Notes: Morphological-type classification from the Marino et al. (20 10) compilation; spectral slope 9 
from {FUV - NUV) color as ^ - [2 + u.4 {t u V '^^UV)/ log(1520/2300)]; 
absolute B maenitude Mb from Bt from RC3 (d^Wucouleui^^l. 1991), m^pt UGC 6022, 
which is from UGC (Nilsor^^TJ)), after correctiotnSnSBHSiSrJ^^ tootnote[T5T and distance, 
the latter derived fronHTJ^ffl?afrtedshift of the parent group, assuming a pure Hubble flow with 
Ho — 75 km s"-*- Mpc"-*-; galaxy stellar mass in Mq from M/ L^, according to Buzzoni 
^200^; SFR in Mq yr"-*- from the dereddened NUV flux, according to footnote |T6l birthrate b 
assumes a current galaxy age of 13 Gyr. 



to be sized up relying on the Galex {FUV - NUV) color. The 
latter directly r elates to the 6 sp ectral slope of galaxy UV emission 
assuming, with lCalzettil l ll999l) . that F{\) oc A*. The intervening 
action of dust works in this sense of absorbing stellar luminosity 
at shorter wavelength making the spectral slope to flatten, so that 
— > 0. O n the other hand, population synthesis models (see e.g. 
iLeitherer et al. 1999; Buzzoni 2002) firmly agree in predicting a 
characteristic spectral slope 6o — —2.1 in case of any dust-free 
star-forming activity driven by the presence of high-mass {M > 
5 Mq) stars. Therefore, the amount of the 6 flattening (i.e. A6 = 
9 — 6o) leads to a straightforward measure of the monochromatic 
attenuation at the different photometric bands 

Internal reddening appears to be quite high in the UV bands, 
with the galaxies loosing typically 2-3 mag at 1500 A, a fig- 
ure which is consistent wi t h the mo re exhaustive studies of 
iRifatto. Longo. & Capaccioli' (1995') and lOordon. Calzetti. & Witll 
1 19971) . Once accounted for reddening and distance, UV abso- 
lute luminosities directly translate into actual SFRol3 while ab- 
solute B magnitude can eventually relate to galaxy stellar mass 



By adopting the ICalzettil il999 ) attenua tion curve, for the magni- 
tude dimming one derives (cf. e.g. Isuzzoiil [20021 A{B) = 1.31 A0, 
A{NUV) = 2.0 Ae and A{FUV) = 2.47 Ad, assuming FU V and NUV 
wave bands to be centered at 1520 A and 2300 A, respectively jNeff et all 
|2008|) . 

1^ According to iBuzzonil j2002l) . the inferred SFR (in Mq yr^^) from 
Galex NUV luminosities is log SFR = log -Fcalcx + 27.70, providing 
to express flux density, in ergs~^ Hz~^. This calibration is for a Salpeter 
IMF with stars between 0.1 and 120 Mq. 



(M*^i ) through the appropriate M/L ratio from the ISuzzonil 
l l2005b template galaxy models, as in Table |2] The galaxy char- 
acteristic birthrate can also be computed from our data, as 6 = 
SFRoto/AI'^i, by assuming to = 13 Gyr. 

The results of our analysis are collected in Table |5] and sum- 
marized in Fig.[T3] lust as a reference, in each panel of the figure we 
also reported the typical error figures of the data (see the ellipses in 
each plot); in particular, we estimate a ±0.2 dex uncertainty on the 
inferred value of log SFRo (mainly as a consequence of the inher- 
ent uncertainty on the value of 9), and a ±0.2 dex uncertainty for 
the derived galaxy mass (accounting for the reddening uncertainty 
and M/L calibration). This eventually leads to a A log b = ±0.3 dex 
for our birthrate estimates. 

From both plots, a nice correlation (p — 0.86) is in place be- 
tween star-formation properties and galaxy (stellar) mass. Indeed, 
the agreement might even be better if one could compare consis- 
tently with the disk stellar mass alone 

A nominal fit to the data in the upper panel of Fig.|13|provides 



log SFRo = 0.56 log A/*^i - 5.6 [Mq y" 
±7 ±7 



(47) 



with (j(logSFRo) = ±0.31 dex. As, by definition, SFRo cx 
6A/gai, the deviation from the one-to-one slope implies a depen- 
dence of the birthrate on the galaxy mass such as 



b(xM. 



gal 



(48) 



as confirmed, indeed, by the point distribution in the lower panel. 
In its essence, the inverse correlation of birthrate with galaxy mass 
closely deals with the nowadays recognized effect of "down-sizing" 
dCowie et alJl996l ; lGavazzi. Pierini. & Bosellill99i) , the imposing 
paradigm in theory of galaxy formation. 

The claimed link with galaxy morphology is another issue to 
be considered, of course. In this case, however, a poorer correla- 
tion (p = 0.59) readily appears for the same data as far as the 
SFRo distribution versus the RC3 classification type "T" is dealt 
with. This is partly in consequence of a less tuned relationship be- 
tween galaxy mass and morphological type, as shown in Fig. II4I 
In addition to the iMarino et alj l l2010l) galaxies, we als o joined 
to the figure the sample of standard spirals studied by iGametll 
(200i), complemented by the low-surface brightness galaxy sam- 
ple of iKuzio de Narav. McGaugh & de Blolirj2004) . and the dl 



collection by Saviane et alj 1 2008 ). For the differen t sets of data 



we de rived M'^i with a similar procedure as for the IMarino et al. 
(201(f) data, namely based on the integrated B (for lGamettll2002l 



and Kuzio de Narav, McGaugh, & de Blok 2004) and H luminos- 
ity (for lSaviane et al.l l2008), respectively. Again, a global fit to the 
data can be adjusted, tentatively in the form: 



T = 10.1 - (A/g*,i/10«) 



(49) 



While the general trend is evident from the data (p — 0.83), 
one could easily verify that any attempt to use eq.(|49t for pre- 
dicting galaxy morphology based on the mass, alone, is a nearly 



While the UV emission mainly traces star formation activity in the 
galaxy disk, our stellar mass derives, on the contrary, from the integrated 
B magnitudes, thus including the bulge photometric contribution. For this 
reason, especially for Sa/b spirals, the value of log M*^^ may be oy eresti- 
mating the disk mass up to a factor of ~ 4 (cf. Table 3 of lBuzzonll2005l) . 
thus leading also to a correspondingly lower value of b. 
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Figure 13. The inferred sta r-formation properties of the spiral-galaxy sam- 
ple by iMaiino et al.' ('2010), referring to three Local-Group Analogs from 
the LGG catalog Garcia ( 1993), as marked and labelled on the plot. Marker 
size indicates galaxy morphology along the RC3 T-class sequence (the 
biggest markers are for Sd/Im galaxies with T ^ 10, while the small- 
est are for Sa systeins with T ~ 1). Current SFR derives from Galex 
N UV fluxe s , dere ddened as discussed in the text, and calibrated according 
to lBuzzonil j2002l) . On the contrary, Galaxy stellar m ass is computed from 
the relevant M / Lb ratio, as from the lBuzzoiiil j2005l) template galaxy mod- 
els. The birthrate assumes for all galaxies a cun'ent age of 1 3 Gyr The small 
ellipses in each panel report the typical uncertainty figures, as discussed in 
the text. All the relevant data for this sample are summarized in Table |5] 
Note the inverse relationship between b and as in the down-sizing 

scenario for galaxy fonnation. 



hopeless taskrn In particular, to a closer analysis, one has to re- 
port that grand-design spirals (type Sbc or class T ~ 4) seem to 
set along a wider spread of mass (nearly two dex in \ogM*^i), a 
feature that partially blurs the otherwise cleaner T -mass relation- 
ship, as already pointed out by [Roberts & HavnesI il994). More- 
over, this picture cannot easily accomodate even dwarf ellipticals, 
that however see m to markedly char acterize t he high-density cos - 
mic environment jCellone & Buzzo ni 2005; G avazzi et al1l2010h . 
As well known, both these features directly call for a long-standing 
quest in the extragalactic debate about the role of primordial gene- 
sis a nd the environment "nurture" to modulate galaxy morphology 
(e.g. 'Dressier" igSC/; 'Thuan, Balkowski. & Tran Thanh Vanl ll992i : 
[eianton et al. 2005; Cooper et al. 2007; Tasca et al. 2009) 

As far as galaxy metallicity is concerned, an int e resting 
pic ture can be devised, relying again on the iGametll ( l2002h 
and lKuzio de Narav. McGaugh. & de BlolJ j2004h contributions. In 
Fig. [15] we further elaborated their original data displaying, in the 
two panels, the mass fraction of fresh gas (including atomic and 
molecular Hydrogen), and the representative [O/H] abundance of 
Hll regions. For the reader's better convenience, a summary of 
these data is also reported in Table |6] The gas fraction C/tot = 
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Figure 14. The observed relationship between morphological RC3 class 
T and galaxy stellar mass according to four dif ferent galaxy sam- 
ples from the work of Gamett (2002) (dots), Mari no et al.r(l2010h (di- 
amonds), Kuzio de Narav. McGaugh. & de Blokl |2004) (triangles) and 
Savianeetal. (20^ (pentagons) (see Table \5\ and \6\ for detai l s). Th e 
value of M* I derives from B ph o tometr y for the Gamett (2002h, 
iKuzio de Narav. McGaugh. & de Blo3 |2004) and ^rino et al. i lOm 
galaxies, and from H absolute magnitudes for the ISaviane et all bom 
sample, throught the appropiiate Al / L ratio, according to lBuzzo rili l2005h . 
Equation (49) provides a fair representation of the data, as displayed by the 
dashed curve. As already pointed out also bv lRoberts & HavnesI i 19941) one 
has to remark, however, the notable dispersion in mass among Sbc galaxies. 

Mga_s/{Mgaa + Afgai) cau casily be derived from the data in force 
of the tight Mgaa/Mgj^i relationship: 



log Afgas = 0.60 log Mg*al + 3.51, 

±8 ±87 



(50) 



with CT(log Afgas) = ±0.25 dex, and p = 0.81. 

The upper panel of the figure shows the somewhat un- 
escapable consequence of the enhanced birthrate among low- 
mass systems. A higher value of b implies that a larger frac- 
tion of fresh gas must still be available at present time. Actu- 
ally, one may even speculate that such a delayed mass process- 
ing is in fact the natural explanation also of the low-metallicity 
figures that characterize, all the way, the dwarf-galaxy population 
in th e Universe ( Arimoto & Tarrab 1990; Lee, Grebel, & Hodge] 
I2OO3I) . This scenario seems to find direct support once exploring, 
as in the upper panel of Fig. [T6] Hll metallicity over th e whole 
range of galaxy mass, for instance by extending again the lGarnet3 
I2OO2) sequence for "standard" spirals with the low-mass samples 
of Kuzio de Naray, McGaugh, & de Blok ( 2004) and Saviane et al] 
( I2OO8I) (see also Table|6j. 

On the other hand, a far more instructive experiment can be 
carried out with the help of eq. ( 131) . Facing a change of the refer- 
ence value for the gas fraction, t/tot, the equation allows us to com- 
pute the induced change in the value of Zgas, by accounting for the 
different ISM dilution factor. Accordingly, the observed gas metal- 
licity could easily be rescaled to any other reference value of Qiai 
starting from the observed value of Gobs- As far as the Hll metallic- 
ity is considered, the expected offset to the observed [O/H] abun- 
dance can therefore be written as 
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The rms of the eq.lED fit can be written as cr(T) ~ (8-0.8 T). MO/H] log |^ 



— log 



l-Gt. 



(51) 
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Table 6. Inferred physical prop e rties for the 'Gamett' 
iKuzio de Narav, McGaugli. & de BloM <2004h . and Saviane et al. 
galaxy samples 



Galaxy ID 


Morph. 
Type ' ' 


log ^^^al 
[Mq] 


ytot 


[0/H]t^^ 






NGC 224 


3.0 


11.44 


0.03 


0.21 


-0.72 


NGC 253 


5.1 


10.80 


0.09 


0.05 


-0.33 


NGC 300 


6.9 


9.95 


0.40 


-0.27 


0.15 


NGC 598 


5.9 


10.22 


0.17 


-0.35 


-0.42 


NGC 628 


5.2 


10.92 


0.25 


0.07 


0.20 


NGC 925 


7.0 


10.52 


0.29 


-0.24 


-0.04 


NGC 1232 


5.0 


11.30 


0.10 


0.05 


-0.29 


NGC 1637 


5.0 


10.21 


0.17 


0.25 


0.16 


NGC 2403 


6.0 


10.31 


0.28 


-0.40 


-0.22 


INLiL Z44z 


3.7 










NGC 2805 


6.9 


10.88 


0.24 


-0.36 


-0.25 


NGC 2903 


4.0 


10.83 


0.09 


0.43 


0.05 


NGC 3031 


2.4 


11.27 


0.02 


-0.01 


-1.03 


NGC 3344 


4.0 


10.30 


0.13 


-0.08 


-0.32 


NGC 3521 


4.0 


10.94 


0.12 


0.29 


0.03 


NGC 362 1 


6.9 


10.67 


0.31 


0.15 


0.40 


NGC 4254 


5.2 


11.17 


0.20 


0.27 


0.26 


NGC 4258 


4.0 


11.18 


0.07 


0.08 


-0.44 


NGC 4303 


4.0 


11. 26 


O.ll 


0.03 


-0.28 


NGC 4321 


4.1 


11.36 


0.13 


0.31 


0.07 


NGC 4395 


8.9 


9.47 


0.51 


-0.57 


0.05 


NGC 5033 


5.1 


11.22 


0.22 


-0.09 


-0.04 


NGC 5055 


4.0 


11.01 


0.14 


0.16 


-0.04 


NGC 5 194 


4.0 


11.20 


0.11 


0.29 


0.00 


NGC 5236 


5.0 


10.77 


0.22 


0.23 


0.29 


NGC 5457 


6.0 


11.12 


0.26 


-0.07 


0.08 


NGC 6384 


3.6 


11.52 


0.10 


-0.13 


-0.51 


NGC 6744 


4.0 


11.45 


0.16 


0.46 


0.36 


NGC 6946 


5.9 


10.97 


0.14 


0.04 


-0.13 


NGC 7331 


3.9 


11.51 


0.08 


0.19 


-0.25 


NGC 7793 


7.4 


9.66 


0.18 


-0.15 


-0.20 




K.UZ10 de 


Narav, McGaugh, & 


de Blok 


(2004) 




F563-1 


8.0 


9.34 


0.49 


-0.81 


-0.23 


F571-5 


9.0 


9.08 


0.77 


-0.75 


0.38 


UGC 1230 


8.7 


9.64 


0.68 


-0.92 


0.01 


UGC 5005 


9.9 


9.27 


0.35 


-0.79 


-0.46 


UGC 9024 


3.0 


10.59 


0.50 


-0.48 


0.12 


UGC 12695 


8.8 


9.81 


0.71 


-0.89 


0.10 


F415-3 


9.9 


8.85 


0.65 


-0.77 


0.!0 


F469-2 


8 7 


9.23 


0.59 


-0 94 


-0 18 


F530-3 


5.0 


10.26 


0.26 


-0.97 


-0.82 


F56I-1 


8.3 


9.50 


0.36 


-0.71 


-0.36 


F563-V1 


9.7 


8.83 


0.45 


-1.26 


-0.75 


F563-V2 


9.7 


9.49 


0.61 


-0.73 


0.07 


F568-6 


7.0 


11.13 


0.23 


-0.19 


-0.11 


F577-V 1 


8.0 


9.59 


0.61 


-0.43 


0.37 


F611-1 


10.0 


8.53 


0.61 


-1.06 


-0.26 


F746-1 


10.0 


10.08 


0.32 


-0.49 


-0.22 


UGC 5709 


7.0 


10.35 


0.40 


-0.03 


0.40 


UGC 6151 


8.8 


9.30 


0.42 


-0.77 


-0.31 






Saviane et al. (2008) 






ESQ 347-G017 


9.0 


8.08 0.67 


-0.91 


0.00 


UGC A442 


8.6 


7JQ 


0.74 


-l.ll 


-0.05 


ESO348-G009 


10.0 


7.62 


0.76 


-0.83 


0.27 


NGC 59 


-3.0 


8.31 


0.62 


-0.76 


0.06 


ESO473^G024 


10.0 


7.06 


0.84 


-1.33 


-0.01 


AM0106-382 


9.0 


6.93 


0.85 


-1.18 


0.19 


NGC 625 


9.3 


8.70 


0.54 


-0.75 


-0.08 


ESO245-G005 


9.9 


8.02 


0.68 


-0.91 


0.03 


DD042 


9.8 


8.08 


0.67 


-1.08 


-0.17 


DD053 


9.9 


7.18 


0.82 


-0.91 


0.36 


UGC4483 


10.0 


6.72 


0.88 


-1.30 


0.15 



(a) : morphological class T, as reported by Hyperleda ^Paturej^^^OOJ); 

(b) : gas fraction Qtot — -^-^gas/ (A/gas + M ) ; mHIie^^^^^aL ^200^ sample Q is 
extrapolated from eq. j50l . 

(c) : for the Sun we assume \og{0/H)Q ^ 8.83 (Grevesse & Sauval 1998) . 

(d) : [O/H] abundance rescaled to Qtot — 0.2. 
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Fi gure 15. Mass fraction of gas and Oxygen abundance of Hi I regions for 
the 



Gametllf2003) (dots) and lKuzio de Narav. McGaugh. & de BlokH2004l) 
(triangles) galaxy samples of Table |6] Marker size is propotional to the 
galaxy morphological class T (the biggest markers are for Sd/Im galax- 
ies with T ~ 10, while the smallest are for Sa systems with T 1). Gas 
fraction accounts for both atomic and molecular phase. Note the inverse re- 
lationships of the data versus galaxy stellar mass, the latter as derived from 
the relevant M / L ratio (see, again, Table|6]for details). 



By entering eq. l l51t with the values of Gtot from Table |6] we 
can eventually rescale galaxy metallicity to the same reference 
gas fraction, say for example C/rof = 0.2, such as [0/H]o,2 — 
[0/H]ohs + A[0/H]. The result of our excercise is shown in the 
lower panel of Fig.[T6]for the joint sample of galaxies. 

Definitely, once accounting for the increasing fraction of fresh 
gas with decreasing galaxy mass, one sees that most of the observed 
trend of [O / H] with Af '^i may simply be recovered by the dilution 
effect of processed mass on the ISM. 



4.4 The Roberts time & the fate of spiral galaxies 

As discussed in Sec. 2.1, the rate at which the mass is chemically 
enriched inside a galaxy follows a sort of universal law; calcula- 
tions show in fact that the burning efficiency factor, F (see Table[Tll, 
is nearly insensitive to the galaxy star formation history. At present, 
about 11% ± 2% of the processed mass in galaxies has been con- 
verted into heavy elements and, according to eq. ([8), the process 
may go ahead for a supplementary factor of (O.ll)"^'^'''^'^ ~ 10* 
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Figure 16. Upper panel: Ox ygen ab undance from Hll regions for the galaxy 
sampl es of iGametll j2002h (dots). iKuzio de Narav, McGaugh. & de Blokl 
i2004h (triangles) and lSaviane et alj J2008h (diamonds). A nice correlation 
is in place with galaxy stellar mass, where low-mass systems appear richer 
in gas (see Fig. 115) and poorer in metals. Once accounting for the differ- 
ent gas fraction among the galaxies, and rescale Oxygen abundance to the 
same gas fraction (Gtot = 0.2, as in the indicative example displayed in 
the lower panel), note that most of the [O / H\ vs. trend is recovered, 

leaving a neai'ly flat [O/ H] distribution along the entire galaxy mass range, 
as predicted by the nearly constant yield-metallicity of the systems. 
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Figure 17. The lifetime of spiral galaxi es as active star-forming systems, 
better recognized as the iRobertsI ( 1 19631) time, is computed along the full 
mass range of the systems. The definite consumption of the primordial 
gas reservoir, as in the original definition of the timescale, Ir, derives 
from eq. 1541 . by parameterizing with the current fraction of returned stel- 
lar mass, /. Present-day galaxies are assumed to be one Hubble time old 
(iff ^ 13 Gyr, see the left scale on the vertical axis). If returned stellar 
mass is also included in the fuel budget, to further extend stai' formation, the 
the "extended" Roberts time, derives from eq. J55t . In both equations, 
the distinctive birthrate relies on the observed down-sizing relationship of 
eq. j47t , while the gas fraction available to present-day galaxies is assumed 
to obey eq. {50\ . eventually exfi'apolated for very-low mass systems (i.e. 
M|^j < 10^ Mq) (dashed curves in the plot). Note, for the latter, that they 
seem on the verge of definitely ceasing their star formation activity, while 
only high-mass galaxies (Af*^j > 10^ Mq) may still have chance to ex- 
tend their active hfe for a supplementary 30-40% lapse of t jj . See text for 
a full discussion. 



of the Hubble time[3 At that time, Hydrogen will definitely van- 
ish even inside stars, and a maximum metallicity ^max has to be 
reached by the system for (1 + i?).^max = 1 — V^. In terms of 
the solar figures, galaxy chemical evolution will eventually end up 
with 

(ymax,^max) (3Y0,1OZq). (52) 

On the other hand, much earlier than the full process comple- 
tion, other outstanding events could mark the chemical history of 
a galaxy. In particular, if one recalls that S (X t^, and takes the 
Hubble time {Ih) as a reference, then the expected evolution of the 
fresh-gas content in the disk scales as 



(1-g) 
(1 - Qh) 




If we set t/ = and recall eq. ( |30t . then the timescale for the galaxy 
to consume its primordial-gas reservoir is 

Curiously enough, the matter annihilation that accompanies nuclear 
burning processes may also result, in the long term, in a not quite negli- 
gible effect for the total mass budget of a galaxy. Previous figures show, 
for instance, that the Mitky Way has already lost in luminosity about one 
thousandth (i.e. a factor of 0.11 X 0.0078, see footnote[2) of its total mass. 
This is roughly the equivalent mass of the whole globular-cluster system 
currently surrounding our galaxy. 



being C/tot and / evaluated at the (cu rrent) Hubble time. After a 
tu time, the so-called "Roberts time' ' ( lRoberts|[T963l) , star forma- 
tion can only proceed further in a galaxy at cost of exploiting its 
own processed mass. Since long, the problem of a reali stic esti- 
mate of the Roberts timescale has been i nvestigated t Sandagel 19861 : 
iKermicutt. Tamblvn. & Congdorj[l994l) . leading to the quite puz- 
zling conclusion that Ir may not substantially depart from one 
Hubble time. As a consequence, present-day galaxies might be 
all on the verge of exhausting their own gas resources. To over- 
come this somewhat embarassing conclusion, it has been empha- 
sized that, after all, even at time Ir a galaxy is not completely gas 
depleted, still counting on a residual buffer A/gas — fM*^i of pro- 
cessed mass (see eq. |31t still rich in Hydrogen and able, in princ i- 
ple, to sustain galaxy SFR for the Gyrs to come dKennicutlll 19981) . 

We can take advantage of our original approach to galaxy 
chemo-photometric evolution to attempt a more confident assess- 
ment of the Roberts time, also in its "extended" definition (tRe), 
such as to account for the "extra-time" provided by processed- 
mass recycling. First of all, according to its standard definition, 
the value of tR can easily be computed with our data along the 
different galaxy mass by means of eq. i50\ . Once getting the 
value of Afgas from the fitting equation, it is immediate to derive 
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Qtot = Mgas/(Afgas + M*^i). Galaxy birthrate follows, as well, 
by relying on eq. ( 147 l l. and recalling that b ~ SFRjf tjf /M'^i. 
The results of eq. ( I54t are displayed in Fig.[T7] maintaining the net 
value of the retumed-mass fraction (/) as a free parameter in our 
calculations. Evidently, current evolutionary scenario can typically 
be sustained by spiral galaxies only for a further 20% of the Hubble 
time. For the most massive systems this is 2-3 Gyr ahead from now. 

When returned stellar mass is included in our budget, too, then 
fresh stars can still form for a supplementary time lapse such as to 
lead to t/tot ~ 0. Recalling again, eq. BOt . then the "extended" 
Roberts time (tfle) follows from eq. ( |53t as 



tRe = 



tH 



tR 



(55) 



again with t/tot and / referring to the boundary conditions of 
present-day galaxies (see Fig.ll7t. Two important conclusions can 
be drawn from the analysis of these results. 

i) For any realistic assumption about the returned mass frac- 
tion, most of the galaxies more massive than ~ 10^^ Mq have 
tR i$ tn- As a consequence, they may already have exceeded their 
Roberts time, being their current SFR mainly driven by processed 
mass rather than by primordial gas. According to the morphology 
vs. mass relationship (see again Fig. [141 this seems to be the case 
for most of Sa/Sb spirals in the present-day Universe. 

ii) When explicitely included in the gas budget, the processed- 
mass contribution certainly adds supplementary chances to high- 
mass spirals for continuing their star formation activity for a few 
Gyrs. However, the relative life extension is critically constrained 
by the balance of the two antagonic processes of mass return 
(/) and mass consumption (b). In absolute terms, eq. ([55} shows 
that the "extended" Roberts time, Ire, is eventually driven by the 
"down-sizing" mechanism, which constrains the b vs. mass rela- 
tionship. The conclusion is that only a marked reduction of the 
present brithrate might allow spirals to sustain star formation far 
beyond one Hubble time. In any case, only a negligible "bonus" 
is awarded to low-mass systems (namely to Sd/Im galaxies less 
massive than ~ 10^^ Mq). For them, in fact, we definitely have 
tRe. ~ ~ in, and just a couple of Gyrs, at most, still seem to be 
left ahead. 



5 SUMMARY & CONCLUSIONS 

In this paper we tried an alternative approach to the study of chemo- 
photometric properties of late-typ e galaxies. F ollowing previous 
theoretical efforts in this sense (e.g. |Pagelll997l) . our analysis relies 
on a basic criterion of energetic self consistency between chemical 
enhancement of galaxy mass, through nuclear processing, and the 
corresponding luminosity evolution of the system. 

This plain physical constraint provides us with a brand new 
interpretative tool, as nearly all the leading theoretical codes along 
the last decades rather approached the problem of galaxy chemi- 
cal evolution in t erms of the s um of individual elemental contri- 
butions (see, e.g. lTinslevlll980l) . This required, in other words, to 
explicitely account for the chemical yields supplied by stars along 
the different range of mass. Clearly, this "analytical" approach is 
the unescapable way when we want to trace the specific enrich- 
ment histoiy of single chemical elements. On the other hand, as 
far as the global metallicity of a stellar aggregate is concerned, the 
analytical process risks to carry along and magnify the combined 
uncertainty from stellar evolution theory making the emerging pic- 



ture to strongly depend on the model assumptions, often through a 
plethora of free tunable parameters. 

On the contraiy, in any "synthetic" approach, galaxy chemi- 
cal evolu tion derives fro m ot her primary physical con straints. The 
works o f ISchmidtl ( [ 19631) and lPagel & Patchetj d 1 9751) are two out- 
standing examples in this sense, assessing the chemical history of 
a stellar system as an implied consequence of its star-formation 
law. As a further variant on this line, our reasoning takes the move 
from the fact that a fixed amount of energy must always be re- 
leased when matter is chemically enhanced through nuclear pro- 
cessing inside stars. Theoretical nucleosynthesis predicts, in fact, 
that 6.98 10^* ergs are produced (see Sec. 2) for each gram of Hy- 
drogen to be converted to heavier elements through non-explosive 
burning processes. A further amount of 1.03 10^* ergs per gram 
should then be added if we want to account for the explosive nu- 
cleosynthesis dealing with the SN events. Overall, it is clear that 
a tight relationship has therefore to be expected between chemical 
enhancement of a galaxy and its past photometric history. 

The theoretical fundamentals which provide the reference 
framework for our anal ysis have b een discusse d in Sec. 2. In par- 
ticular, we rely on the iBuzzonj ( Il989l Il995h original code for 
stellar population synthesis, further elaborated to build up three- 
zone galaxy template s along the Hubble morphological squence 
( lBuzzonill2002[ . |2005|) . For its special relevance within the global 
properties of late-type systems, our study is mainly focussed on the 
disk evolution, which broadly marks the macroscopic look of spiral 
and irregular galaxies. The standard scenario of the present calcu- 
lations considers a Salpeter IMF, with stars formed within a mass 
range between 0.1 and 120 M©. In addition, metals and Helium 
proceed from processed stellar mass assuming a fixed enrichment 
ratio i? = 3, as discussed in Sec . 2.3. 

Following iBuzzoriH ( |2005[) . an important assumption of our 
models is that stellar birthrate is a distinctive property of galaxy 
morphology and mass (the latter, in force of the observed T-M*^i 
relationship, as shown in Fig. II4t . Consi stently with the down- 
sizing scenario for galaxy formation (iCowie et al.ll9%h . b is there- 
fore s een to decrease, in our models, with increasing M*^i dCavazzil 
Il993h . As a consequence, SFR scales with time according to a 
power law, and it is not directly coupled to the mean gas density 
of the disk as im posed, on the contrary, by other classical schemes 
( ISchmidtll95^ . 

Chemical enhancement inside a galaxy is effectively assessed 
by the burning efficiency factor, as displayed in eq. ([§} of Sec. 2.1. 
Table [T] reports a notable property of F, being it nearly insensitive 
to the galaxy star formation history. Along a full range of possible 
evolutionary scenarios, in fact, we always end up at present time 
with F — 0.11 ± 0.02. This means, in other words that, within 
just small (±20%) individual differences, galaxies performed very 
similarly everywhere in the Universe as efficient engines for ener- 
getic exploitation of the barionic matter through nuclear synthesis. 
Accordingly, at any epoch we can univocally define a reference fig- 
ure, that we called "yield metallicity" (Z, see eq. [24|in Sec. 2.3), 
which traces the mean composition of the stellar mass. 

The role of SNe in the more general context of disk chemical 
evolution has been assessed in some detail in Sec. 2.2. A compar- 
ison of the energetic budget involved in quiescent and explosive 
nucleosynthesis indicates that SNe are extremely powerful mass 
processors; for fixed amount of released energy, in fact, they bum 
a factor of seven more mass than "normal" stars. Following yet 
classical arguments (e.g. lMatteucci & Greggioiri986l) we know that 
the component of high-mass stars (M > 5 Mq) enables a very 
quick chemical enrichment during early evolution of galaxies. This 
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"prompt" chemical enhancement is the result of the explosion of 
both core-collapsed objects, in the form of type II SNe, and ac- 
creted WDs as a result of binary interactions, like in type la SNe. 

Yield metallicity is therefore abruptly raised to log Z ~ —3 
on a timescale of a few 10^ yrs, a figure that further increases by 
~ 0.15 dex within the first Gyr of galaxy life (see Fig. |2). As, by 
definition, the cumulative action of "prompt" SNell and la must be 
proportional to the galaxy stellar mass i\/gai, the net effect of their 
contribution is to set a steady offset (that we quantified in Sec. 3.2 
in [Fe/H] ~ —1.10 dex) to galaxy yield metallicity. 

In this framework, for a Salpeter IMF, the implied yield 
metallicity evolves as 2 oc (see eq. |25). For the previ- 

ous arguments, this can be seen as a nearly universal law, only 
marginally dependent on the galaxy birthrate (see, for instance. 
Fig. |3). An immediate consequence of this important property 
is that IMF, rather than SFR, is the key player to set metallic- 
ity at primeval epochs, when t ^ and SNe were governin g 
chemical evolution (e.g. lScannapieco, Schneider. & Ferrarall2003h . 
In particular, non-Salpeter stellar populations, strongly biased to- 
ward high-mass stars as in a flatter IMF, have to be invoked to fit 
with the peculiar presence of bright super-solar galaxies at high 
redshift, as sometimes observed (e.g. .Maiolino et al., ,2006t) and 
also envisaged by some updated the oretical schemes (see, fo r in- 
stance, the interesting explorations o f lKroupa & Weidneil200l and 
iKoppen. Weidner, & Kro"uD3l2007l) . Figure [T] is an illustrative ex- 
ample of this mechanism when tracing the contribution of SNell to 
metal enhancement with varying the IMF power-law index. 

A major advantage of taking yield metallicity as a reference 
marker in our analysis is that Z can easily be related to other key 
physical processes that constrain galaxy chemical evolution in its 
different facets. The ISM metal enrichment is an evident example 
in this sense. Chemical composition of the gaseous phase inside 
a galaxy is the resulting balance of two basic mechanisms. From 
one hand, one should carefully consider the bulk of processed mass 
returned by stars to the ISM through SN events and quiescent stellar 
wind. On the other hand, this input has to be properly assessed in 
the physical context of dynamical evolution of the disk, such as 
to trace the way processed mass is eventually "diluted" within the 
fresh primordial gas. 

In tackling the problem (Sec. 3.1), we tried to maintain our 
formal treatment as much as possible free from any "ab initio" 
assumptions that might bias our conclusions. The returned mass 
flow is therefore dealt with in eq. ( 128) in terms of the parameterized 
fraction f{t) of net stellar mass returned to the ISM within the 
time t (see footnote [TOt. In our notation, this component is always 
traced separately from the primordial gas component G{t), so that 
the global gas amount in the disk comes throughout as a sum of the 
two nominal contributions, as in eq. ( I30t . In its formal elaboration 
(see eq. |29|and|31|), the ISM metallicity (^gas) can eventually be 
set in terms of yield metallicity, and derives from just two leading 
parameters, namely Q and / (see Fig.|7}. For the latter, a safe upper 
limit has been placed, such as / < 0.3 at any age for a Salpeter 
IMF, according to a full discussion of the SSP case, as in Sec. 3.1.3. 

The case of the Milky Way must be the obvious and preem- 
inent testbed for our theoretical framework. In particular, the ex- 
pected evolution of the yield metallicity is compared, in Fig. [4] 
with the observed AMR, as traced by different stellar samples in 
the solar neighborhood. By definition, at every epoch Z caps the 
value of the ISM metallicity, and therefrom edges the composition 
of newly formed stars. As discussed in Sec. 3.1, the study of the 
A[Fe/H] ~ log(Zgas/-Z) difference is richful of important infor- 
mation to constrain the Galaxy star formation history. We have been 



able, for instance, to set a firm upper limit to the Galaxy birthrate, 
such as b < 0.5 (see ea.l33t. while the current contribution of the 
primordial gas cannot exceed roughly one-third of the local mass 
density of the disk (see eq. [35). A further boundary condition can 
also be posed to the chemical enrichment ratio, that must be i? < 5. 

A similar comparison, performed between "analytical" model 
predictions from several reference codes in the literature and ob- 
served AMR, is summarized in Fig.|5] The broad range of envis- 
aged scenarios predicted by theory leads to a spread of up to 0.5 dex 
in the value of stellar [Fe/H] at a given age. We want to look at 
this figure as a measure of the intrinsic limit of theory to master 
the combined interplay among the many leading mechanisms that 
modulate chemical evolution. To a closer analysis, however, it has 
to be recognized a clear tendency of models to predict a sharper 
chemical evolution of Galaxy disk, such as to exceed, in most cases, 
the observed stellar metallicity at the present epoch. Evidently, the 
shallower enrichment rate indicated by the observations points to 
a still unsettled problem with the primeval steps of disk evolution, 
and more specifically with the pre-enrichment processes related to 
the well known G-dwarf problem. It is interesting to see, in this 
regard, that when the observed stellar [Fe/H] distribution is ex- 
plicitely dealt with in our analysis, as in Sec. 3.2.1, the implied 
birthrate (b = 0.6 ± 0.1) tends to oversize the corresponding figure 
derived from the AMR. Such a more "delayed" star formation is in 
fact a sign for primeval stars to have formed "elsewhere" (alias in 
the "thin" disk, as claimed by the standard established scenario). 

An outstanding dyscrasia of any theoretical scheme for Galaxy 
chemical evolution is that no spread in [Fe/H] can be admitted 
among coeval stars populating the same region of the system, be- 
ing the AMR univocally settled by the boundary physical condi- 
tions of the model. This is evidently at odds with what we observe, 
at least in the solar neighborhood, where an intrinsic a[Fe/H] up 
to 0.2 dex seems an unquestionable feature for the local stellar pop- 
ulation (see Fig.[8j. As an important consequence of the metal dis- 
persion, one may actually conclude that, disregarding age, about 
four out of five stars in the solar vicinity approach the expected 
yield metallicity within a factor of two. When interpreted in terms 
of mixing properties between enriched and primordial gas, such a 
tuned distribution clearly demonstrates that star formation in the 
Galaxy only proceeded, all the time, in a highly contaminated en- 
vironment. This claim can be better quantified by means of Fig.|7] 
which confirms that the processed stellar mass returned to the ISM 
must in fact be the prevailing component to gas density in star- 
forming regions if we want [Fe/H] of fresh stars at any epoch to 
closely match (i.e. within A[Fe/H] ~ 0.3 dex) the yield metallic- 
ity, as observed. 

This evidence severely tackles the classical scheme a la 
Schmidt, where star-formation strength is fully driven by just the 
mean gas density. On the contrary, what it might be the real case, 
is that star formation process in spiral galaxies behaves like a wild 
"outbreak", that stem from previously contaminated regions to pro- 
gressively spread across the rest of galaxy body (lNepveulll98"l. 
1989; Barnes 2004). 

The possible implication of the Milky Way scenario to the 
more general case of late-type galaxy evolution along the Hubble 
morphological sequence has been the focus of Sec. 4. As far as 
closeby external galaxies in the local Universe are concerned, ob- 
servations still hamper any firm assessment of galaxy metal abun- 
dance. As a matter of fact (see Fig.lI2t, the galaxy-to-galaxy scat- 
ter of [Fe/H] distribution seems to largely overcome the nominal 
accuracy of any individual estimate, a feature that makes evident 
the intrinsic difficulty in defining a fair "representative" metallicity 
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when such a composit e stellar environment is dea lt with, as in the 
late-type galaxy disks jArimoto & JabIoiikalll99lh . 

On the other hand, the combined study of the other leading 
parameters can add important pieces of information, that greatly 
help constraining the problem. In this regard, our discussion along 
Sec. 4.3 basically grounds on three relevant relationships, that 
emerge from the observational properties of different galaxy sam- 
ples. In particular, 

i) the down-sizing mechanism appears to govern star forma- 
tion in the local Universe, as well, with a clear relationship between 
current SFR and galaxy stellar mass (see Fig. 113): 

ii) the "delayed" star formation among low-mass galaxies, as 
implied by the inverse h-M*^i dependendence, naturally leads to a 
more copious gas fraction when moving from giant to dwarf galax- 
ies, as observed indeed (see Fig.llSt: 

Hi) the physical relationship of star formation mechanisms 
with galaxy stellar mass can also aproximately fit into a straight 
morphological scheme (Fig. II4t . where galaxies with decreasing 
mass more likely take the look of later-type spirals. It has to remain 
clear, however, that mass not morphology is the primary pa rameter 
that govern late-type galaxy properties ( Boselli et al.booi]) . 

Once accounting for galaxy chemical properties, by taking for 
instance Oxygen abundance of Hll regions as a tracer, the blurring 
trend with morphology leaves space to a much cleaner relationship 
with galaxy stellar mass, as in the upper panel of Fig. [16] Facing 
the fact that low-mass systems also appear richer in gas and poorer 
in metals, one might even speculate on a simple scenario where the 
observed trend is in fact merely the result of the softening mecha- 
nism, that leads enriched stellar mass to mix into fresh unprocessed 
gas in the ISM. Actually, when galaxy data are properly rescaled to 
a fixed reference value of t/tot, through eq. ( 15 It , one sees (lower 
panel of Fig. [T6l ( that most of the [O / H] vs. M*^^ trend is re- 
covered, leaving a nearly flat [O/H] distribution along the entire 
galaxy mass range. This result is consistent with our theoretical 
predictions for a nearly constant yield metallicity along the entire 
late-type galaxy sequence. 

A striking evidence, when sampling the local Universe, is that 
spiral galaxies appear all dangerously close to shut down their 
venture as stars-forming systems, being on the verge of exhaust- 
ing their own gas resources. A so special case, and our somewhat 
suspicious temporal location as privileged observers, di d not es- 
cape t he analysis of many authors in the literature, leading [Roberts! 
( Il963l^ first to wonder about the real "residual life" left to spiral 
systems according to their current gas capabilities. Relying on our 
notation, the so-called Roberts time {tn) easily derives in its orig- 
inal definition, namely as the timescale for the primordial gas to 
vanish (see eq. [54}. On the other hand, and more interestingly, 
one would like to assess the even more general case, when the 
global gas consumption (i.e. including returned stellar mass) even- 
tually sets (or overcome) the ultimate lim it to galaxy star formation 
( iKennicutt, Tamblvn. & Congdon|[T994h . 

Equation J55t shows that, as far a the down-sizing prescrip- 
tions set the reference relationship between birthrate and galaxy 
mass, one has to conclude that spiral galaxies cannot escape an age 
where no stars will form at all. In fact, if returned mass is left as 
the only gas supplier to the ISM, then the maximum birthrate that 
a galaxy can sustain is of the order of femax — f tn / {1 ~ f) ^ 
//(I — /), which must evidently be much less than ~ 0.45, for a 
Salpeter IMF. As a result, only massive 

(A^gai > 10" Mq) Sa/Sb 
spirals may have some chance to extend their active life for a while 
(namely for ~ 0.3 tn, as indicated by Fig.|I7|(. On the contrary, no 
way out seems to be envisaged for dwarf systems, that will soon 



cease their star formation activity unless to drastically reduce their 
apparent birthrate below the 6max figure. 
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